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ABSTRACT 

Trees provide environmental, economic, and social advantages in urban areas. Knowing 
the extent and location of tree canopy in a municipality is an important step in quantifying these 
benefits. Spatial and temporal tree canopy analysis was performed for the city of Columbus, 
Georgia, by categorizing the National Agriculture Imagery Program (NAIP) aerial imagery for 
2005, 2010 and 2015 into tree versus non-tree land cover type using unsupervised classification 
procedures. Air pollution removal rates from the i-Tree program were applied to this evaluation 
providing an estimate of the city’s tree air quality benefits. The city’s canopy overall has 
remained steady at 52% of the 38,143 hectares that compose the municipality for the years 2005 
(89% accuracy), 2010 (93% accuracy) and 2015 (93% accuracy). Percent tree canopy within the 
city’s 53 census tracts ranged from 13 to 75%. Tree loss due to development in south central, 
north, and north-eastern areas was offset by forest regrowth, having been cleared prior to 2005. 
These trees remove 1,700 tonnes of five critical air pollutants (CO, NO2, O03, PM2.5, PMio, and 
SOz) and sequester 256,000 tonnes of COz2 annually, based on i1-Tree’s first-order valuations. 

Since trees influence fine particulate matter (PM2:s) and the health impacts of PMz.s are 
great, a second study was conducted to better understand how tree stand formation controls 
PM2.s. Three portable, fine PM sensors (AirBeams) were used among three tree canopy 
configurations (dense tree buffer, n=5; small tree line, n=6; and U-shaped, n=4) to determine if 
stand design effects PM2.s concentrations in open areas near trees. AirBeams were evaluated and 
found to have reliability, ease of use, repeatability among units, and stability across the study 
period. Overall results between open and tree concentrations were not significantly different. 
Site by site observations indicated that dense tree buffers (3 of the 5 sites) trap PMz2.s5 resulting in 


higher tree particulate concentrations in the buffer zone and small tree lines (5 sites) had no 





effect on PM2.s. U-shaped tree stands interactions are dependent on location of the open area 
within the tree stand in relation to notable PM sources. While wind direction played a role in 
particulates reaching sampling locations, proximity to and type of PM source had the largest 
impact on local PM2.s concentrations. 

Urban canopy cover recommendations are made so cities can benefit from ecosystem 
services that trees provide, but simply adding trees does not mean these benefits are fully 
utilized. Tree type, tree design, and tree placement, i.e. in available space and proximity to 
pollution source, need to be considered. Utilizing high spectral imagery and low-cost, portable 
sensors can help cities determine the best tree placement and design to aide in air pollution 


reduction. 
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INTRODUCTION 

Trees provide many benefits in urban areas, such as improving air quality, sequestering 
carbon dioxide (Nowak & Crane, 2002), filtering water (Booth, 2005), and decreasing urban heat 
islands (Bolund & Hunhammar, 1999). These environmental benefits also have health 
advantages, like improving senior longevity (Takano, Nakamura, & Watanabe, 2002), lowering 
the number of autism cases (Wu & Jackson, 2017), and lowering mortality with cleaner air 
(Tiwary, 2009). People value trees mostly for their shade, air quality improvements, and 
“calming effects” (Lohr, Pearson-Mims, Tarnai, & Dillman, 2004). The economic gains to a city 
can include everything from lower crime rates (Kuo & Sullivan, 2001) to increased housing 
prices that generate property tax revenues (Anderson & Cordell, 1988; Donovan & Butry, 2010). 

Researchers have taken different approaches to understanding the values of trees. Top- 
down approaches involve quantifying trees regionally by classifying high spectral imagery into 
landcover classes including tree canopy cover (Nowak, 2012). Canopy is defined as the area of 
land covered By tree leaves, trunks, and branches as seen from an aerial perspective (Northern 
Research Station, 2017). Once canopy is analyzed, models can be applied to the findings to 
estimate tree benefits (McPherson & Simpson, 2002; Nowak & Crane, 2000). Conversely, field 
studies aimed at quantifying tree benefits in a locale are preferred, especially when few field 
studies exist to support models that apply air quality values to trees (Pataki et al., 2011; Setala, 
Viippola, Rantalainen, Pennanen, & Yli-Pelkonen, 2013). 

While all ecosystem services and economic attributes provided by trees are important, 
this research focuses on the removal of air pollutants by trees. The first investigation looks at the 


value of applying broad level air quality assessment to cities using high spectral canopy analysis 


and air pollutant removal rates, obtained from the i-Tree Tool (www.i-Tree.org), for Columbus, 





Georgia, a municipality in western Georgia, USA. The air pollutant removal rates using the i- 
Tree Tool are based on first-order approximations, with environmental and meteorological data 
from one location often accounting for several counties in a region. As such, local research 
should assist in determining services and possible disservices of tree placement or removal 
practices (Nowak & Greenfield, 2008; Nowak et al., 2014). While other limitations exist in 
using the i-Tree model, fine particulate matter (PM2:5) is the one pollutant that poses additional 
specific limitations. The i-Tree model’s PMz.s uses a positive removal rate in counties with low 
wind and a negative removal rate (meaning increase in PM2:s) in counties with high wind and 
low rain (Hirabayashi, 2014). Trees are a temporary resting location for PM2.5, and local weather 
conditions, especially wind and precipitation, can resuspend particles into the air or bring to the 
ground (Nowak et al., 2014). As a result, it was determined that the second part of this research 
would focus on tree and fine particulate matter interactions by tree buffer arrangement, as this 
interaction is more complex than other air pollutants removed from the air column by trees. 

Particulate matter (PM) is among the six criteria air pollutants regulated by the 
Environmental Protection Agency (EPA) under the National Ambient Air Quality Standards 
(NAAQS) as part of the Clean Air Act (Girard, 2014). PM2.s5 causes major health related issues 
when inhaled (Hemon & Fechner, 2014), and has been linked to over 100,000 deaths annually in 
the United States (Fann et al., 2012). 
Study Goals and Scope 

As the city of Columbus was the first city in the state to become Tree City USA certified 
(and has remained so for 39 years; GFC, 2012), maintaining a working knowledge of the tree 


canopy and its services (or disservices) in the community is vital. The goal of this research is to 


pair spatial and temporal analysis of canopy with air quality monitoring to quantify tree benefits 





and aide in future tree planning for the city of Columbus. The final product will contain: a 
thematic tree canopy map and percentage breakdown of tree canopy by census tract for 2005, 
2010 and 2015; tree canopy change over the ten-year period; a first-order estimation of 
Columbus air quality benefits; and the results of a field study examining tree effects on local 


PM2.s5 levels within tree canopy stands and adjacent open areas. 





CHAPTER 1 “i SPATIAL AND TEMPORAL CANOPY COVER ANALYSIS 
1.1 Introduction 

Municipalities for decades have focused on vegetative planning to improve water 
filtration, reduce air pollutants, support economic growth, and abate climate change (Howard, 
1965; Miller, 1988; Platt, Rowntree, & Muick, 1994; Young, 2010; Escobedo, Kroeger, & 
Wagner, 2011; Roy, Byrne, & Pickering, 2012). Urban planning often incorporates these 
ecosystem services into local urban design with vegetation in mind, but determining what to 
account for, whether air pollution abatement, social improvements, or economic values take 
priority, is complex (Thomas & Geller, 2013). Knowing the amount and location of tree canopy 
in an urban environment is a mandatory first step as municipalities plan for future development. 

The term canopy, for the purposes of this research, means the area of land covered by tree 
leaves, trunks, and branches as seen from an aerial perspective (Northern Research Station, 
2017). US cities and counties produce Urban Tree Canopy (UTC) assessments using high 
spectral aerial or satellite imagery to classify land cover thereby ascertaining tree canopy amount 
and distribution. The USDA Forest Service’s Northern Research Station and the University of 
Vermont’s Spatial Analysis Laboratory created procedures and have assisted cities in developing 
UTC assessments to help with urban tree planning (Northern Research Station, 2017). 
Municipalities use similar procedures by classifying high spectral imagery to assess tree canopy 
spatially and temporarily in order to enhance “green” planning. As an example, Atlanta, 
Georgia, completed a UTC assessment using satellite imagery (2-foot pan-sharpened, 4-band 
data) through the Georgia Tech Center for Geographic Information Systems in 2014. That 
assessment determined the city had 47.9 percent canopy cover in 2008 (Giarrusso & Smith, 


2014), 





Land cover classification using high spectral imagery is a top-down approach to 
determine tree canopy versus the bottom-up approach used when surveying individual trees in an 
area. Each approach has its advantages and disadvantages, with the top-down approach being 
the best approach for assessing amount and location of trees in larger areas the size of 
municipalities (Nowak, 2012). The use of imagery to assess different land cover types is based 
on the idea that dissimilar objects, like water, vegetation, and roads, have unique spectral 
signatures because they reflect and absorb wavelengths of electromagnetic radiation (EM) 
differently (Kachhwaha, 1983; Keranen & Kolvoord, 2014). For example, water absorbs red and 
near-infrared wavelengths (0.76—0.90 um), while vegetation reflects these wavelengths (0.63-— 
0.90 um). This difference allows the two land types to be distinguished using multispectral 
imagery and image analysis programs, like ESRI ArcGIS (Fox, 2015). Multispectral imagery is 
composed of three to seven bands of pixels with values 0 to 255, lower values are darker and 
higher values are lighter. Each band represents either visible or infrared wavelength ranges, i.e. 
for Landsat images band 1 is visible blue (0.45 to 0.52 um), band 2 is visible green (0.52 to 0.60 
Lm), band 3 is visible red (0.63 to 0.69 um), band 4 is near infrared (0.76 to 0.90 um), band 5 is 
short-wave infrared | (1.55 to 1.75 um), band 6 is thermal infrared (10.4 to 12.5 um), and band 7 
is short-wave infrared 2 (2.08 to 2.35 um). These bands are combined and analyzed based on 
known spectral signatures of objects to classify an area of concern (Keranen & Kolvoord, 2014; 
Fox, 2015). 

Once a municipality’s tree canopy is known, ecosystem services (i.e. air quality, water 
filtration, and reduction of urban heat) can be estimated using tree models, like i-Tree Tools 


(online tools and software developed by US Forest Service, Davey Tree Expert Company, 


National Arbor Day Foundation, Society of Municipal Arborists, International Society of 





Arboriculture, and Casey Trees) and ArcGIS based CITYGreen (software developed by 
American Forests with rates based on Urban Forest Effects (UFORE) methods, an earlier version 
of i-Tree). These models estimate UTC effects on air pollutant removal and carbon dioxide 
sequestration and storage. These removal rates were established through modeling the 
combination of tree canopy across the United States, leaf area index values, pollution removal 
rates by trees given local pollutant concentrations, and pollutant deposition rates based on local 
meteorological data. For i-Tree, the monetary value of these ecosystem services was applied 
based on health incidences and associated costs that would be avoided with pollutant removal 
(Nowak, Hirabayashi, Bodine, & Greenfield, 2014). The UFORE model took a similar 
approach, but used fewer cities and applied monetary values based on prevented health and 
tourism loss (Nowak & Crane, 2000). 

The i-Tree Tool are a good starting point for tree planning and quantifying associated 
benefits of city trees. The i-Tree Tool incorporates data from across the United States and 
assesses removal rates for five of the six criteria air pollutants (carbon monoxide - CO, nitrogen 
dioxide - NO2, ozone O3, particulate matter — broken into PM2.s and PMio, and sulfur dioxide - 
SO2). The model also applies sequestration and storage rates for carbon dioxide (CO2). 
Empirical studies quantifying tree impacts on air quality are limited (Pataki et al., 2011), and 
models projecting tree reductions of air pollutants range from 0.13 percent for PM2.s removal in 
urban settings to 0.51 percent for ozone removal in rural settings (Nowak et al., 2014). Other 
research has quantified the ability of trees to reduce air pollutants, but the focus is on only one or 
two criteria air pollutants (reviewed in Nowak et al., 2014). The air pollutant removal rates 
using the i-Tree Tool are based on first-order approximations, with environmental and 


meteorological data from one location often accounting for several counties in a region. One 





Italian study found good agreement when 1-Tree ozone removal rates were compared to local 
level field measurements (Morani et al., 2014). 

The Natural Resources Spatial Analysis Laboratory conducted a Georgia Land Use 
Trends analysis including tree cover for the state of Georgia. Landsat satellite data (30 m x 30 m 
resolution) was used to create GIS databases for the state with an overall accuracy of 85 percent 
for years 1974, 1985, 1991, 1998, 2001, 2005, and 2008 (Kramer, 2016). The study found an 
eight percent tree canopy loss in Muscogee County between 1991 and 2005 (GFC, 2012). The 
city of Columbus (consolidated with Muscogee County, Georgia) did not have a recent 
assessment of tree canopy nor a high resolution assessment needed for city wide tree planning. 

Maintaining a working knowledge of the Columbus city tree canopy and its services (or 
disservices) in the community is important given the potential environmental implications. This 
research seeks to fill the knowledge gap regarding the city’s tree canopy through spatial and 
temporal analysis of its UTC. The goal of this study is to assess the city of Columbus tree 
canopy and estimate its associated air quality benefits. The knowledge gained will aide in 
providing sound recommendations to the city on advantageous locations for future tree planting 
and removal and will enhance planning and policies concerning development with vegetation in 
mind. 

Disparity in tree canopy can be found in a look at the spectral imagery of Columbus, 
Georgia. It is visibly apparent that a greater, healthier canopy exists in the northern portion of 
the area. The National Agriculture Imagery Program (NAIP) imagery is a rich green color above 
highway 80 (denoted in red in Figure 1). Conversely, tree canopy is scarce in the downtown 
industrial area (the portion of the city that appears greyish-white in the mid-west portion of the 


image). This contrast in percent canopy is consistent with the American Forests 





recommendations for ideally 15 percent canopy in downtown and industrial and 50 percent 
canopy in suburban residential areas (American Forests, 2002). When disparity in canopy exists 
across the city, the environmental benefits of the urban trees are also unequal city-wide. Tree 
benefits are applied at a small scale given the spread of atmospheric conditions and sources of air 
pollutants in urban settings (Nowak & Greenfield, 2008; Tyrvainen, Pauleit, Seeland, & de 

Vries, 2005). 

This research intends to answer the following question: will aerial imagery analysis, 
quantifying tree canopy spatially and temporally, highlight large tree canopy and air quality 
benefit disparities over time across the city of Columbus? It is hypothesized that a disparity will 
exist between tree canopy within the study domain, with northern areas of the municipality 
having the most canopy and downtown having least canopy, resulting in disproportionate air 
quality benefits across the city. Additionally it is hypothesized that the tree canopy will decrease 


over the time period examined. 
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Figure 1. City of Columbus service region (excluding Fort Benning) 2015 1-meter, 4-band NAIP 
natural color image with inset map showing location of Columbus, Georiga. 








1.2 Methods 

1.2.1 Study Area — Columbus is located along the western border of Georgia, USA 
(Figure 1 inset, 32° 29’ 32” N, 84° 56’ 25” W). The city and county (Muscogee) governments 
are consolidated, therefore, the area of the city is the County land area of 56,045 ha (138,490 
acres). Part of the Fort Benning Army base is located in southeastern Muscogee County. 
Excluding this portion of the County (17,902 ha), the city of Columbus has a land area of 38,143 
ha (94,253 acres). The landscape to the north-northeast is dominated with agriculture and pine 
forest found throughout the southeastern United States, while the south-southwestern landscape 
is urban. Columbus is the second most populous municipality in Georgia with population of 
189,885 in 2010 (U.S. Census Bureau, 2010). 

1.2.2 High-Resolution Imagery - The 2010 and 2015 tree canopy were analyzed using 
the 1-meter ground sample distance (GSD) spatial resolution, 4-band National Agricultural 
Imagery Program (NAIP; 1m x Im spatial resolution) imagery of Muscogee County. The 2005 
tree canopy was analyzed using 2-meter GSD spatial resolution, 3-band NAIP imagery, as this is 
what was available. NAIP produces digital orthoimages roughly biannually (Georgia imagery 
exists for years 2001-2002, 2005-2007, 2009-2010, 2013, 2015, and 2017; USDA, 2017) by 
aerially photographing agricultural regions during the growing season, usually between July and 
September. The 1-meter GSD spatial resolution, available for free through the USDA program, 
offers the best resolution publicly available for Columbus in recurrent years, and it is, therefore, 
the best available imagery of the city accessible to city planners and other researchers for future 
tree canopy analysis. The 2005 NAIP imagery was obtained through the Columbus City 
Planning Department in compressed county mosaic format (Figure 2A). The 2010 city NAIP 


imagery was obtained from the U.S. Department of Agriculture Farm Service Agency, Aerial 
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Photography Field Office in Digital Ortho Quarter Quad (DOQQ) tiles format containing 25 
separate DOQQ files (Figure 2B). The 2015 NAIP imagery was available online to download 
via the Aerial Photography Field Office in compressed county mosaic format (Figure 1). All 


three were projected in the UTM coordinate system, North American Datum of 1983. 
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Figure 2. City of Columbus natural-color image A) 2005 2-meter spatial resolution, 3-band 
NAIP and B) 2010 1-meter spatial resolution, 4-band NAIP. 


Imagery classification was conducted using ESRI ArcGIS 10 software. NAD 1983 State 
Plane Georgia West FIPS 1002 Feet was used as the projected coordinate system as requested by 
the city of Columbus GIS Division. The city service boundary shapefile was used to clip Fort 
Benning from the imagery. Imagery was gathered at different times of day and on different days. 


The 2010 imagery was flown between two separate months. As a result, shadows, which 
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complicate classification, are at different angles in different sections. Recognizing flight pattern 
reduces this potential classification error. Quarter quadsacre closely followed the flight pattern 
used to gain the imagery. For this reason, the clipped city image was split using the DOQQ 
shapefiles to reduce classification error. The 2010 imagery was provided in 25 DOQQ tiles, so 
analyzing each separately was the most effective analysis technique (Figure 3A). For the 2015 
imagery, DOQQ shapefiles were combined into four large images to group based on similar 
topographical sections (e.g. urban versus forest) and to reduce processing time (Figure 3B). The 
2005 imagery, as received, was not cleanly mosaiced with north-south alignment issues across 
the image (Figure 3C highlights this issue). To reduce error, this imagery was cut along these 


mosaiced sections and analyzed in 10 sections (Figure 3D). 
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Figure 3. A) 2005 NAIP Columbus natural-color image showing an area of misalignment. The 
sections used to analyze B) 2005, C) 2010, and D) 2015 NAIP imagery. 
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1.2.3 Tree Canopy Image Classification - For larger areas that an analyst is not fully 
acquainted with, performing unsupervised classification can reduce error (Rozenstein & Karnieli, 
2011). Therefore, for each section of the NAIP imagery, an unsupervised iso cluster 
classification process was conducted clustering the image bands (3-band for 2005, 4-bands for 
2010 and 2015) into 40 classes. The 40 classes were visually interpreted and assigned labels of 
tree or non-tree (see Appendix A for iso cluster values by year). Postprocessing procedures 
involved mosaicing the classified sections into one raster of the whole city. The Majority Filter 
and Boundary Clean tools were used to clean the image by filling in areas of no data and 
smoothing the edges of tree and non-tree clusters (Keranen & Kolvoord, 2014). This approach 
yielded three classified thematic maps spatially showing tree canopy across Columbus, Georgia, 
in 2005, 2010, and 2015. 

1.2.4 Classification Accuracy Analysis — A simple random sampling scheme with a 
sample size of 500 reference points for each year (i.e. 1,500 total points) was used to assess 
accuracy. Using multinomial probability theory, the following equation was used to determine 


sample size: 


_, Brj(1-7j) 
ar ee 
bj 


N 
where B is the upper a/k percentile of chi square distribution with | degree of freedom, k is 
number of classes (2), a is acceptable error (0.05), 7t; is the proportion of trees in the 
classification (0.52), and b; is confidence interval and precision (0.05). The multinomial model 
and simple random sampling satisfy assumptions of the kappa statistic (K-hat), which is used to 
calculate the significance of the error matrix table generated during the accuracy assessment 


(Congalton, 1991). Using this method, 1,500 reference points (500 for each of the three thematic 


classification rasters being assessed) were randomly generated using the Create Radom Points 
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tool in ArcGIS (Figure 4). All three years were assessed at each point with a separate value 
applied for each year based on landcover type in that particular year. Assessing all 1,500 points 
for all three years increases sample size, which in turn increases confidence in the classifications 


(Dicks & Lo, 1990). 
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Figure 4. Position of 1,500 randomly generated reference points used to check accuracy of 
classified thematic maps with in the Columbus, Georgia, service region. 





The reference imagery data used to access classification accuracy for 2005, 2010, and 
2015 was Google Earth™ imagery. This tool was selected over actual ground truthing at the 
physical location because landcover changes rapidly. Google Earth™ serves as a good reference 
tool as the imagery within the tool has high-resolution and offers the historical images needed for 
assessment (Gongalton, 1991; Olofsson et al., 2014). The history bar within Google Earth was 
utilized to access imagery from 1993 to 2017 for the region. The tool allows for rotating views, 
viewing imagery from different angles, and street view, which helped in determining tree versus 
non-tree when shadows were prominent. All years of imagery available (1993, 2003, 2005, 


2006, 2007, 2009, 2010, 2011, 2012, 2014, 2016, 2017) were used to assess the accuracy of the 
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classification, especially the years before and after those of interest (1.e. 2005, 2010, 2015), as 
described in Olofsson et al. (2014). 

All 1,500 reference points were manually assessed using Google Earth™ and given a 
value of 1 for tree and 2 for non-tree. The Extract Values to Points tool was used in ArcGIS to 
compare reference points to the classification raster for each of the three years (see Appendix 1, 
Table 4). Accuracy was then assessed using the error matrix table and corresponding Kappa 
coefficient (Congalton, 1991). 

1.2.5 Spatial and Temporal Analysis — The Spatial Analysis tools in ArcGIS (Price, 
2014) were used to determine percent tree canopy of the 53 census tracts using the 2015 US 
Census TIGER tract shapefile for the city of Columbus service region. A comparison of tree 
canopy change over time by tract was conducted and assessed by evaluating percent change by 
census tract between the three classified thematic maps. Changes within tracts were further 
evaluated to determine cause of any differences found, i.e. tree loss due to development and 
timber harvesting or gains due to tree plantings. 

1.2.6 Air Quality Benefit Analysis - The i-Tree Tool was used in conjunction with the 
tree classification results to estimate urban tree canopy tree canopy air quality benefits. This tool 
applies average air pollutant removal rates and monetary values based on county level data. 
These county level rates were determined by combining tree canopy analysis, leaf area index 
(LAI) values, pollution removal rates by trees given local pollutant concentrations, and pollutant 
deposition rates based on local meteorological data (Nowak et al., 2014). The 2001 National 
Land Cover Database was used to determine tree cover and percent of cover that was evergreen, 
while the LAI values were found using the MODIS/Terra global Leaf Area Index product. Tree 


removal of air pollutants was determined using a statistical model that combined total tree cover, 








evergreen percentage, LAI, local weather, and local air pollutant concentration data 
(Hirabayashi, 2014). Monetary value was estimated based on health incidences and associated 
costs that would be avoided with pollutant removal (Nowak et al., 2014). Table 1 contains the 
removal rates derived using this process for Columbus, Georgia. 


Table 1. Tree air pollution annual removal rates and related monetary values for Columbus, 
Georgia, using 1-Tree, developed by USDA Forest Service (Nowak et al., 2014). 


b 
(Removed annually) (tonnes/hectare-year) ($/tonnes) 






$40.25 
$39.00 


$39.00 


1.3 Results 

1.3.1 Classification Accuracy Assessment — A 93 percent overall accuracy was found for 
the 2010 and 2015 classifications, while the 2005 classification had an accuracy of 89 percent 
(Table 2). The user and producer accuracy are the same for both tree and non-tree for the 2010 
thematic map, so error was spread evenly between error of omission and commission (Table 2). 
The user accuracy (error of commission) for trees, i.e. the percent of trees correctly classified, is 
highest in the 2015 classification at 95 percent, with 2005 also being good at 92 percent. The 
kappa statistics for all three years is relatively high showing good agreement between reference 
data and thematic map data after accounting for agreement by chance. User accuracy for non- 
tree is 4 percent lower than tree for 2015 and 6 percent lower for 2005 classifications. The 


producer accuracy (error of omission) for trees, i.e. the percent of pixels correctly labelled as 


trees, is 3 percent lower for tree versus non-tree for the 2015 classification and 4 percent lower 
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for 2005. In summary, the 2015 and 2005 classifications represent actual referenced trees better 
than non-trees, while the percent of pixels correctly labelled as non-tree is higher. 


Table 2. Error matrices for 2005, 2010, and 2015 classifications containing user and producer 
accuracy, k statistics (78, 86, and 87 percent respectively), and overall accuracy (89, 93, and 93 
percent respectively) results. 


2015 Accuracy Assessment 2010 Accuracy Assessment 2005 Accuracy Assessment 
Reference Reference Reference 
Data Data Data 
Thematic Tree Non- | Map User's Tree Non- | Map User's Tree Non- | Map User's 
Map Data Tree | Total | Accuracy Tree | Total | Accuracy Tree | Total | Accuracy 
l 










Non-Tree aon | 713 659 | 712 103 | 619 | 722 


7 
Total 
Accuracy 
Overall Accuracy = 93% Overall Accuracy = 93% Overall Accuracy = 89% 
K-hat = 87% K-hat = 86% K-hat = 78% 


1.3.2 Spatial and Temporal Dynamics of Tree Canopy Coverage — In 2015 and 2010, 
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the city of Columbus tree canopy covered 52 percent of the area, equivalent to 19,815 and 
19,809 ha (48,964 and 48,949 acres), respectively. In 2005, the tree canopy made up 53 percent 
of landcover, equivalent to 20,012 ha (49,453 acres). 

In 2015 the tree canopy cover in the 53 census tracts ranged from 13 to 75 percent of land 
cover. The range was 10 to 75 percent in 2010 and 9 to 73 percent in 2005 (Figures 5 and 6). 
While the overall canopy coverage for Columbus remained steady (2005-2015), the change over 
time within certain tracts and in certain areas of the city is notable (Figure 7). In 2005, 5 tracts 
had less than 20 percent canopy. This number dropped to 2 tracks with less than 20 percent 
canopy in 2015. The number of tracks with 20 to 39 percent canopy changed from 18 in 2005 to 


26 in 2015. The tracts within the 40 to 59 percent canopy range decreased from 24 in 2005 to 21 


in 2015. The tracks with the highest canopy (60 percent and over) decreased from 6 in 2005 to 4 
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in 2015 (Figure 6, see Appendix A for Table 5 summarizing by tract tree canopy and Table 6 by 


tract air quality benefits). 
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Figure 5. City of Columbus thematic tree canopy map for A) 2005, B) 2010, and C) 2015. 
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Figure 6. Percent tree canopy by census tract: A) 2005 ranging from 9 to 73 percent UTC, B) 
2010 ranging from 10 to 75 percent UTC, and C) 2015 ranging from 13 to 75 percent UTC. 
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Figure 7. City of Columbus tree canopy change by census tract between 2005 and 2015. Light 
green represents losses (22 tracts) and dark green represents gains (13 tracts) in canopy over the 
ten-year period. 


1.3.3 Air Quality Benefit Analysis - Approximately, $4 million in health-related savings 
can be attributed to the removal of 1,700 tonnes (1,900 tons) of pollutants annually by trees 
(Table 3). $10 million of savings is due to 256,000 tonnes (282,000 tons) of carbon dioxide 
sequestered annually by Columbus trees (Table 4). Additionally, the 20,000 hectares (49,000 
acres) of trees store 5.6 million tonnes (6.2 million tons) of carbon dioxide valued at $218 
million (i.e. this is a long-term value). 


Table 3. Columbus, Georgia, tree air quality benefits using USDA Forest Service i-Tree tool 
(Nowak et al., 2014). 
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Table 4. Columbus, Georgia, tree carbon dioxide sequestration and storage using USDA Forest 
Service i-Tree tool (Nowak et al., 2014). 
2015 2010 2005 


Polutanks Columbus a Columbus ae Columbus | 2005 Monetary 
(Removed 


Monetary Monetary 
Removal Removal Removal Value ($/yr) 
annually) (tonnes/yr) Value (/y?) (tonnes/yr) Value @iyr) (tonnes/yr) 


256,300 | $9,995,634 | 256,221 | $9,992,572 | 258,860 | $10,095,460 
5,583,416 | $217,751,913 | 5,581,706 | $217,685,205 | 5,639,177 | $219,926,585 


Air quality benefits across Columbus are best visualized by applying the removal rates to 





trees within each census tract (Figures 8 and 9). Trees in the northern portion of the city (tracks 
101.07, 108.02, 102.03, 102.01, and 103.01) remove the largest tonnage of air pollutants per unit 
area. The trees in the downtown areas (the southwestern portion of the city) remove the least, 
with the midtown trees removing slightly higher amounts of pollutants. This trend matches the 
tree canopy across Columbus. Air pollution removal rates through the i-Tree Tool are calculated 


based on tree coverage. 
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Figure 8. Annual air pollution removal by census tract. Numbers (except 11 and 25) represent 
tracts with largest air pollution removal. Numbers 111 and 25 represent tracts with lowest 
pollution removal. 
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Figure 9. Annual CO? sequestration by census tract ranging from 1.7 to 9.7 tonnes/hectare. Dark 
blue represents tracts with largest and light blue the least sequestration. 


1.4 Discussion 

1.4.1 Classification Accuracy Assessment — The lower accuracy seen in the 2005 data is 
due to the NAIP imagery having a spatial resolution of 2-meters and only 3-bands. Using 4-band 
imagery allows for greater distinction between water and vegetation, both of which can have a 
greenish hue especially when water has high nutrient content. Water is not reflective, but rather 
absorbs EM radiation in the near-infrared (band 4) while healthy vegetation is very reflective in 
this band (Fox, 2015). Additionally, the 2005 iso cluster rasters had lower resolution with larger 
clusters covering multiple landcover types, i.e. trees, buildings, and road were in one clustered 
pixel group. Given these constraints, the 2005 NAIP imagery was more difficult to analyze, 
which increased error. 


The main source of error in the 2010 and 2015 classification rasters was due to 


performing an iso cluster unsupervised classification. Distinguishing shadows between those 
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concealing trees and those concealing non-tree landcover was difficult as these were combined in 
at least one iso cluster class for most clipped sections. Often iso cluster classes combined tree 
and non-tree features. For example, tree and non-tree vegetation share at least one or two iso 
cluster classes in each clipped section because spectral signatures for grass, bushes, and trees can 
overlap. Reclassifying these classes as only tree or non-tree increased error. Unfortunately, all 
classification processes have error whether computer generated, as with unsupervised 
classification, or human error, as seen with the supervised classification process due to lack of 
familiarity with the region being analyzed (Rozenstein & Karnieli, 2011). The accuracy of the 
2010 and 2015 classifications is good compared with other classifications of NAIP imagery 
found in literature (Davies et al., 2010; Li et al., 2014; Moskal, Styers, & Halabisky, 2011). Also 
of note, the 2010 NAIP imagery had added error with clouds covering a field of trees in the 
northeast portion of the image. 

1.4.2 Spatial and Temporal Dynamics of Tree Canopy Coverage — The ideal canopy 
cover for an urban area as stated by the American Forests Urban Forest Program is 40 to 60 
percent for forested states like Georgia (Leahy, 2017). The city of Columbus falls well within 
this range when the city is considered as a whole, possessing 52 percent canopy cover between 
2005 to 2015. The 200-hectare difference between 2010, 2015 and 2005 is negligible when 
classification error is considered. While the aggregate canopy cover did not change over that 
time period, the canopy within the 53 census tracts did change over time. These results highlight 
the city’s recent development and forestry practices. 

The greatest loss in canopy over the ten-year period occurred in the lower middle census 


tracts (Figure 7). Based on interpretations of classification changes over time and changes seen 


during the error check, there are two main reasons for tree loss: development and removal of 
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residential trees. The majority of the UTC loss seen over time was associated with development. 
In many cases, trees were cleared between 2005 and 2010, but buildings and pavement were not 
in place until after 2010. These areas of development have lower canopy cover to start with, and, 
therefore, the tree loss leaves a greater impression than areas with greater percent canopy. The 
increase in impervious surfaces in these areas affects to air quality too, as roads and businesses 
increase vehicle traffic to these areas. Ornamental trees are often planted at new businesses and 
shopping areas, but these trees are smaller than the mature trees removed during construction. 
Ornamental trees, like crape myrtles (Lagerstroemia indica L.), have small leaf area indexes and 
mature tree heights, which makes them poorly suited for reducing air pollutants (McPherson, 
Simpson, Peper, & Xiao, 1999; Yang, Chang, & Yan, 2015). 

Census tract 25, the area between the Chattahoochee river and south of Highway 280 
containing the Columbus Civic Center (Figure 8), had the second smallest tree canopy for all 
three years. This tract experienced a 9 percent gain in tree canopy between 2005 and 2015 
because the city arborist and local tree organizations focused on tree plantings in this area (S. 
Jones, personal communication, October 27, 2017). Based on a thorough examination of the 
area and discussions with the city arborist, little can be done to further improve tree canopy in 
this area unless businesses get involved, even with a canopy cover of only 18 percent in 2015. 
Much of the land is owned by the city in the form of public parks with ball fields and parking 
lots. The remainder of the tract is private property. 

Other tree canopy gains over the ten-year period are due in large part to tree growth on 
forestry lands previously cleared for timber. This growth was evident in the northeast corner of 


Columbus (Figure 5). Another noteworthy area is the northeast tract, known as Midland (tract 


101.07). While this area has experienced little overall change in canopy between 2005 and 2015, 
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a lot has changed in canopy location across this large census tract. The most eastern portion of 
the tract has many pine tree farms that were harvested around 2005 and have since been 
replanted. The western to middle portion of this tract has experienced a lot of development with 
the expansion of neighborhoods and businesses at the expense of tree canopy. While the forest 
regrowth offset the losses due to development, another timber harvest would reduce tree canopy 
for this area and the city as a whole. Due in part to the city tree ordinance enacted in 2002, 
Columbus maintained marginal loss of trees despite large gains in impervious surfaces. Key to 
the reduction of tree loss is the mandate that requires new business developments plant trees in 
parking lots. 

When adding trees to an urban environment, planning often focuses on location and types 
of trees to best provide the ecosystem services trees offer. Attention should be given to utilizing 
as many of tree benefits as possible in addition to air pollution reduction, like water 
management, social and recreational values, and noise reduction (Miller, Hauer, & Werner, 
2015; Grey, 1996; Jim, 2004). City owned property, i.e. parks, city buildings, monuments, 
cemeteries, and right of ways, lacking tree coverage is the first priority for planting locations 
(Grey, 1996). Columbus has done a good job of managing trees in many of these areas, but 
downtown municipal buildings lack appreciable tree canopy. In the Columbus downtown area, 
cemeteries, the medical center, and businesses comprise the land available for planting trees. 
Increasing trees in these vegetation sparse areas will involve educating businesses on the value of 
trees. 

When compared with 5 other counties and their associated major cities in the 


Southeastern United States, Muscogee County has the best canopy cover (Table 5). However, 


the urban portion of Columbus only has more canopy cover than Montgomery, Alabama. As 
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Chatham is located on Georgia’s coast, it’s land cover is 32 percent water, which reduces land 
available for tree plantings (Plan-It Geo, 2015). The five counties viewed for comparison, 
except for Chatham, Georgia, have the same trend as Muscogee County: they all possess greater 
canopy cover than their major cities. This alludes to the idea that air pollutants are most likely 
being produced in areas with lower numbers of trees to reduce the pollution. Charlotte, North 
Carolina, has similar population density over land area as Columbus, but Charlotte has a greater 
canopy cover as compared to the urban portion of Columbus. This suggests that Columbus can 


improve its tree canopy in the developed portions of the city. 
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1.4.3 Air Quality Benefit Analysis - Spatially, the city tree canopy differs greatly from 
north to south. The census tracts (101.07, 102.03, 102.01, and 103.03) north of highway 80 and 
tract 108.02 (second most eastern tract below 101.07) are considered the northern portion of the 
city. These four and tract 108.02 (previously Fort Benning land) are not as developed as the rest 
of the city, containing mainly forest and agriculture landcover. These five tracts comprise 53 


percent of the area for the city of Columbus. Only 16 percent of the Columbus population 
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resides in this northern portion of the city. The remaining 48 tracts make up the other 47 percent 
of the land in the southern portion of the municipality. This distinction, north versus south, 
separates the mainly urban, developed portion of Columbus (south) from the agricultural, rural 
portion (north). 

The north portion of Columbus contains two-thirds the city tree canopy. Not surprisingly 
these 5 tracts experience the most air quality benefits of trees (1,141 tonnes of air pollutant 
removal, 170,000 tonnes CO2 sequestered annually). If these 5 northern census tracks were 
removed from Columbus (leaving the urbanized portion of the city), it would only have 37 
percent tree canopy capable of removing an estimated 576 tonnes of air pollutants and 
sequestering 86,000 tonnes of CO2 annually. Many of the city’s shopping centers and businesses 
exist in the south central portion of the city, so most residents must travel within the southern 
section. Therefore, the majority of the air pollutants are being produced (via vehicles and 
businesses) in the portion of the city with the least number of trees. 

Urban forest management involves diversifying types of trees planted and selecting trees 
that can remain stable, improve air quality, and not emit high amounts of volatile organic carbon 
(VOC). VOCs contribute to air pollution and can lead to higher particulate concentrations 
(Miller et al., 2015). The best tree species to reduce pollution are typically unpopular trees to use 
in street and residential planting (Yang et al., 2015; Simpson & McPherson, 2011; Curtis et al., 
2014; Benjamin, Sudol, Bloch, & Winer, 1996). Conversely, popular trees, like oaks, offer great 
air pollutant reduction but also emit high VOCs during spring and summer seasons (Curtis et al., 
2014; Bolund & Hunhammar, 1999). This points to the need to diversify the types of trees 


planted across cities, prioritizing the species with the best overall performance. 
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The methods used here to quantify tree canopy coverage and its associated air quality 
benefits have limitations. Comparing city canopy cover using aerial imagery at five-year 
intervals may not provide adequate time to detect differences at a city-wide scale. Analysis 
every ten years, is better for tracking these changes. However, technological advancements are 
expected during a time period of ten-years, which makes comparable imagery difficult. NAIP 
imagery improved over the 10 years used in this analysis, from 3-band, 2 m resolution in 2005 to 
4-band, 1 m resolution in 2010. Starting in 2017, three states had NAIP imagery available with 
50 cm resolution. It is likely that this higher resolution imagery will be available for all states 
soon (USDA, 2017). Satellite imagery is also improving, offering better resolution and more 
bands than NAIP (WorldView-2: 0.5m resolution with 8 bands). However, these satellite images 
are not free, like NAIP. 

Meneguzzo, Liknes, and Nelson (2013) found the unsupervised approach to NAIP 
imagery classification overestimates tree clusters as compared with object based image analysis 
(OBIA) and better reflects photo-interpreted results compared to ground-based or bottom-up 
approaches. Ground surveying trees in a city on a block by block approach is often the next step 
after quantifying canopy using high spectral imagery (Miller et al., 2015). Tree surveys are 
helpful in identifying tree health, height, and type, which better assists in planning at a street 
level. 

Since the i-Tree Tool is a first-order assessment of air pollutant removal associated with 
trees, it is not possible to estimate the degree of accuracy and variability associated with its 
predictions. The developers of i-Tree acknowledge that there are limitations in using this 
method (Nowak et al., 2014). Removal rates are calculated in part based on air pollution 


concentrations and meteorological data, which are gathered at county and regional levels. For 
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example, SO2, PM2.s, O3, and meteorological data were retrieved from data collected at the 
Columbus Airport (AQS Site ID: 13-215-0008), which is centrally located in Muscogee County. 
PMio data was collected at Cusseta Road Elementary school (AQS Site ID: 13-215-0011) in 
south Columbus. NO2 and CO data were collected near Atlanta, Georgia, (AQS Site ID: 13-089- 
0002) approximately 145 km away from Columbus (EPA, 2015). Thus, the removal rates for 
SO2, PM2.5, PMio and Os are better applied at the city scale, given the data used to determine 
these rates was obtained within Columbus, as compared to the removal rates for NO2 and CO. 
Vehicles and energy production are main sources of NO2 and CO emissions in urban settings 
(Girard, 2014). As Atlanta, Georgia, is more heavily populated with more vehicles than 
Columbus, Georgia, relying on NO2 and CO data gathered in the Atlanta metropolitan area to 
create Columbus removal rates raises concerns. With fewer vehicles in Columbus, it is 
conceivable the pollution concentrations in the area are lower than the Atlanta area. Therefore, 
removal rates for these two air pollutants by trees in the Columbus area are likely overestimated. 
In this study, pollutant removal rates were applied at the census tract level. While useful for 
highlighting variation in the removal of air pollutants across the city, this approach may not 
accurately represent pollutant attenuation. As air pollution is generated locally, the trees in the 
northern portion of the city may not be removing air pollutants at the rates estimated if the air 
pollution does not exist in these areas. 

In addition to the limitations caused by the lack of air pollutant concentration data, too 
few studies quantify tree reductions of air pollutants (Pataki et al., 2011; Setala et al., 2013). 
This scarcity of data is a noteworthy drawback to the i-Tree model. Additionally, the PM2.5 
concentrations removed by trees may be insignificant compared to the levels in the air (Witlow, 


2009). 
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Over half the monetary savings of the five criterion pollutants found using i-Tree are 
attributed to PM2.s removal (4 percent of the 1,700 tonnes of air pollutants removed annually). 
PM25 health concerns have been well researched (Sarnat, Schwartz, & Suh, 2001; Pope III et al., 
2002; Schlesinger et al., 2006; Fann et al., 2012), but the tie between trees interactions with 
PM2.5 and the credited health benefits has not been well researched (Pataki et al., 2011). 
Estimates from the i-Tree model suggest that Columbus trees offer $4 million in health savings. 
This estimate may not be an accurate assessment. 

As an example, consider PM2.s and its complex interactions with trees. This relationship 
cannot fully be captured in a simple removal rate, even if the rate was generated using Muscogee 
County specific data. Other air pollutants such as O23, NOz, SO2, and CO (being gaseous) have 
simpler relationships with trees because they are removed from the atmosphere by leaf stomata. 
In contrast, temperature, wind direction and wind speed alter PMz.5 deposition and resuspension 
(Hemond & Fechner, 2014; Nowak et al., 2014; Tong, Whitlow, MacRae, Landers, & Harada, 
2015). County PM2.s removal rates vary in the i-Tree model and are positive or negative 
(increases in PM2.s) depending on county wind and precipitation conditions (Hirabayashi, 2014). 
While PM2.s rates used in the model come from county specific data, the interaction with trees 
and PM2.s is clearly complicated. This complex interaction between PM2.5 and local tree canopy 
conditions requires additional research. 

While the approach employed in this research to quantify trees’ air quality services has 
limitations, the benefit of this top-down approach in locating and determining change in canopy 
over time allows for city planners to develop tree plans that can improve local environmental 


conditions. Conducting a thorough tree benefit analysis at the street level would be ideal, but, 


access, time, and funding are large hinderances to this endeavor. A standard model, like i-Tree, 
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provides city planners with information that can be used to gain support and funding for future 
tree plantings. As long as caution is taken by weighing the approximated air quality reductions 


and monetary values in light of the model’s restrictions, this approach is helpful to start 


conversations relating to vegetation planning at the city level. 
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CHAPTER 2 — TREE ARRANGEMENT PARTICULATE MATTER TESTS 
2.1 Introduction 

Air pollution reduction benefits by trees in urban areas have been linked to health 
improvements (Beckett, Freer-Smith, & Taylor, 2000; Tiwary, 2009; Nowak et al., 2014). The 
air quality benefits of trees for gaseous air pollutants is a simpler relationship compared to 
particulate pollutants because trees remove gaseous pollutants from the atmosphere through leaf 
stomata. In contract, trees are a temporary resting location for fine particulate matter, PM2:s. 
Tree interactions with particulates is dependent on weather conditions such as temperature, 
relative humidity, wind direction and wind speed all of which alter PM2.5 deposition and 
resuspension (Hemond & Fechner, 2014; Nowak et al., 2014; Tong, Whitlow, MacRae, Landers, 
& Harada, 2015; Cai et al., 2017). While studies have been conducted to better define this 
relationship, studies often discuss the need for testing across regional conditions (Ortolani & 
Vitale, 2016; Nowak & Greenfield, 2008; Nowak et al., 2014; Witlow, 2009). The complex 
interactions between trees and particulates warrants additional examination. 

2.1.1 Particulate Matter and Trees - Particulate matter (PM) is among the six criteria air 
pollutants monitored by state agencies and regulated by the U.S. Environmental Protection 
Agency (EPA) under the National Ambient Air Quality Standards (NAAQS) as part of the Clean 
Air Act (Girard, 2014). PM is categorized by aerodynamic diameter into coarse, PM10 — 2.5 to 
10 um in diameter, and fine particles, PM2.5 - 2.5 um or smaller in diameter. The residence 
time of PM10 is minutes to hours, being removed from the air due to gravitational settling. As a 
result, PM10 travels less than 100 km. PM2.5 has a residence time of days to weeks. It is often 


removed through dry deposition and rain and travels 100s to 1000s of kilometers (Wilson & 


Spengler, 1996). In interactions with trees, PM2.5 levels on leaves are lower than PM10 due to 
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gravitation deposition properties (Beckett, Freer-Smith, & Taylor, 2000; Freer-Smith, 2005; 
Seb¢ et al., 2012). Fine particulates are solid/liquid aerosols created through anthropogenic 
activities e.g., fossil fuel combustion, wildfires, steel making, and natural processes e.g., sea 
spray, pollen, vegetation releasing volatile organic compounds (Hemond & Fechner, 2014; 
Girard, 2014). These microscopic particles are most distressing for health reasons as inhaling 
these aerosols can cause respiratory and cardiovascular complications (Sarnat, Schwartz, & Suh, 
2001; Pope III et al., 2002; Schlesinger et al., 2006). Inhalation of PM2.5 caused an estimated 
130,000 deaths in the U.S. in 2005. For comparison, 4,700 deaths were attributed to ozone 
exposure in the same year (Fann et al., 2012). 

Field studies monitoring PM2.s within and near tree buffers report varying results. 
Several small-scale experiments in the New York City area indicate that tree buffers limit PM2.5 
dispersion causing concentrations to be elevated in close proximity downwind from tree lines, 
while PM2.s concentrations quickly decrease in open areas (Tong et al., 2015). The New York 
City study and two other studies conducted in and near Beijing found PM2.5 concentrations are 
higher within dense tree canopy buffer or forests as compared with open areas (Tong et al., 2015; 
Liu, Yu, & Zhang, 2015; Chen et al., 2015). A Detroit, Michigan field study found vegetation 
barriers caused particulates to decrease more gradually beyond the tree stand than open areas 
(Brantley, Hagler, Deshmukh, & Baldauf, 2014). Witlow (2013) found that particulates 
increased with distance behind tree stands. Another study reported that tree buffers decreased 
particulate matter beyond tree stands when the wind is from the direction of the road (Baldauf et 
al., 2008). These reported results were conducted across different weather conditions in various 
locations, which further points to the need for site specific analysis of tree fine particulate air 


quality benefits. 
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The EPA has a guide for designing vegetation barriers along roadways in order to reduce 
air pollution and recommends denser tree lines to reduce air flow and stop pollutants near the 
street (Baldauf, 2016). In recent years, studies have been conducted in various tree 
configurations. Most studies, discussed above, indicate higher PM2:s concentrations exist within 
denser canopies (Tong et al. 2015; Whitlow, 2013; Tong, Chen et al. 2015). One study modelled 
the impact of air pollutant reduction among six tree designs finding dense tree buffers most 
effective (Baldauf, Isakov, Deshmukh, & Zhang, 2016). A few studies have modelled the tunnel 
effect created by trees lining either side of the street, which traps pollutants between the tree 
lines increasing concentrations (Gromke, 2011; Cai et al., 2017). These studies investigated the 
trees near roadsides, concentrating on the reduction of vehicle produced pollutants. Vehicles are 
a main particulate pollution sources in cities, but other sources like restaurants, prescribed burns, 
and utility companies can also contribute to PM2.s concentrations (Zheng et al., 2002). 

Additional field studies are needed to better understand the relationships between urban tree 
stand arrangements and PM2.s5 concentrations, especially accounting for local PMz2.5 sources and 
atmospheric conditions. 

2.1.2 PM2.5 Instrumentation - Fine particulate matter as regulated by the EPA is 
monitored by states using in situ continuous monitors that are accurate and expensive equipment, 
e.g, the Tapered Element Oscillating Microbalance (TEOM) and b-attenuation monitoring 
(BAM) analyzers (EPA, 2013). These monitors are located in larger cities across the U.S. (EPA, 
2015), and are designed to measure regional PM levels. Air pollution across cities is 
heterogeneous, depending on localized interactions of PM2.s, weather, and vegetation, and small 


portable sensors can be beneficial in identifying areas of focus and concern. Low-cost, portable 


PM sensors can be used to bridge the gap in knowledge between regional particulate data and 
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neighborhood level exposure to pollutants. However, caution should be taken when using these 
portable devices due to accuracy and reliability concerns (Jiao et al., 2016; Lewis & Edwards, 
016; Rai et al., 2017; Snyder et al., 2013; Wang et al., 2015; Manikonda, Zikova, Hopke, & 
Ferro, 2016). 

The EPA encourages cities to study air quality at several locations using a variety of 
sensors to assess local air quality conditions (EPA — Smart City Challenge, Green Cities project). 
With advances in sensor technology, people can monitor their local air quality by operating a 
personal, portable PM device. While not promoting specific equipment, the EPA has supported 
these efforts by providing agency research on portable PM devices and data regarding how 
individuals and communities can use these devices to facilitate urban planning (EPA — Air 
Sensor Toolbox). 

The Community Air Sensor Network (CAIRSENSE) project was conducted by the EPA 
and Georgia Environmental Protection Department to assess the accuracy of a few of these low- 
cost particulate matter devices including the AirBeam. The AirBeam is a low-cost portable fine 
particulate matter sensor developed in 2013. In CAIRSENSE, three AirBeams were tested for 
168 days and compared to federal equivalent method (FEM) monitors. The results indicated 
mid-level agreement (r = 0.65-0.66; Jiao et al., 2016). This result was comparable to another 
low-cost PM device, the Dylos (r = 0.63-0.67). Notably, CAIRSENSE found a strong 
association among the three AirBeam devices tested (e.g. r = 0.99; J ijao et al., 2016). Another 
recent study compared AirBeam performance to reference instruments in field tests and found 
similar results (AirBeams to each other: r = 0.99, to GRIMM 11-R r’ = 0.66 to 0.71; Mukherjee, 


Stanton, Graham, & Roberts, 2017). While the AirBeam has mid-level agreement to FEM, these 


devices have utility as portable field devices when conducting comparative studies. 
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2.1.3 Study Goals - This research project investigated the question of the relative 
relationship of PM2.s concentrations within varying tree canopy types versus adjacent open areas. 
The AirBeams provide a useful tool to address this question. Since the complex relationship 
between trees and particulate matter is dependent on local factors (e.g., particulate sources, 
weather conditions), this research presents an opportunity to test the efficacy of AirBeams in the 
field at different locations and under different atmospheric conditions. Thus, one question posed 
by this research is whether a portable, low-cost monitoring device (specifically AirBeams) can 
be used to effectively compare PM2.5 among differing tree canopy configurations? Based on 
previous AirBeam studies (Jiao et al., 2016; Mukherjee et al., 2017), it is hypothesized that the 
devices will be sufficient for measuring relative relationships among tree canopy configurations, 
but they will not be accurate for quantifying actual PM2.5 concentrations in the field. As such, 
the first goal of this research is to assess the field capabilities of the AirBeam. 

The PM2.5 concentrations within open areas will be compared to adjacent tree buffers 
composed in the following tree arrangements: dense tree buffer (width > 45 m), small tree line 
(width < 30 m), and U-shaped tree arrangement (Figure 10). The following question will be 
addressed in this approach: Do open areas near PM sources differ in PMz2.s concentrations as 
compared with adjacent tree stands of various configurations? It is hypothesized the dense tree 
stands will have higher PMz2.s concentrations than adjacent open areas, small tree line PM2.5 
concentrations will have no appreciable difference to adjacent open areas concentrations, and the 
U-shaped tree arrangements will have highest relative concentrations of PM2.5 within the open 


areas adjacent to the trees. Columbus has average low wind speeds, PM2.5 should be trapped by 


dense trees and disperse less in open areas along streets. Therefore, the primary goal of this 
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study is to examine the relationship between PM2:s and trees by investigating varying tree canopy 


configurations. 
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pire 10. Pichi of neielibotine open areas and an example of A) dense tree buffer, B) small 
tree line, and C) U-shaped tree arrangement with examples of open (blue arrows) and tree 
(orange arrows) sample locations. (Images taken in 2017 via Google Earth™.) 


agle Earth = 


2.2 Methods 

2.2.1 Study City Atmospheric Conditions and PM Sources - Columbus, Georgia, is 
located along the western border of the state (32° 29’ 32” N, 84° 56’ 25” W). The Fort Benning 
Army base is located to the southeast of the city. This city has an annual precipitation of 46 
inches (NWS, 2017). Maximum monthly temperatures range from 57 °F to 92 °F and monthly 
minimum temperatures from 36 °F to 73 °F (lows in January and highs in July and August). 
Monthly relative humidity stays close to the annual mean of 65 percent, peaking at 71 percent in 
August with low of 62 percent in February and March (NCEI, 2017a). The city has mean winds 
speeds of 2.5 m/s with highest winds, 2.9 m/s, in the winter months. Mean wind direction is 
variable throughout the year, with the highest frequency of winds from the east (Figure 11; 
NCE], 2017a). Climatological data used in this study was retrieved from the weather station 
located at the Columbus Metropolitan Airport (WBAN: 93842, Lat/Long 32.5161°, -84.9422°). 


The city’s high annual precipitation and low annual wind speed should be conducive for PMz.s5 to 
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be trapped by trees and brought to the ground rather than resuspended into the air (Hirabayashi, 


2014). 
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Figure 11. Columbus wind rose (data from 2007 to 2016) indicating highest frequency of winds 
from the east. Wind direction frequency is marked by the 2 and 4 percent circles. 


Fort Benning is located to the southeast of Columbus. Many of the city’s worst air 


quality days occur during controlled wildfire burns at Fort Benning. These burns increase PM2:s 


concentrations, among other pollutants. Several studies have monitored the air pollution at 


increasing distances around Fort Benning (Achtemeier, 2011; Baumann, 2005; Liu, 2010; 


Odman, 2012) because this pollution impacts the region. Like other cities, vehicles, residential 


wood burning, and meat cooking are main sources of fine particulate matter (Reff et al., 2009; 


Zheng, Cass, Schauer, & Edgerton, 2002), which have more localized influence as compared to 


the prescribed burns. 
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2.2.2 Instrumentation — In order to quantify PM2.5, this research utilized the AirBeam 
(Figure 12A). This device employs a Shinyei PPD60PV particle sensor and Bluetooth to 
transmit data through a smart phone app called AirCasting (Android app, Figure 12B) or a 
website (aircasting.org/map). The AirBeam can record data while mobile or in fixed position. 
The AirBeam reports PM2.s(ug-m™), temperature (°F), percent relative humidity, and sound 
level (decibels) every second, minute, or hour in real time (Heimbinder & Besser, 2014). The 
Shinyei PPD60PV particle sensor uses the light scattering method to count particulates that cross 
the path of the encased infrared light (Figure 13). This particle sensor has a concentration 
measurement range of 0.5 to 300 ug-m* (Shinyei Technology Co., Ltd). The AirBeam has an 
output resolution of 0.0001 ug-m* when recording at 1-minute intervals and 0.01 wg-m”? when 
recording at 1-second intervals (Heimbinder & Besser, 2014). The data were sent from the app 
to email in comma-delimited format and converted to an Excel spreadsheet for analysis. The 
AirBeams each have a unique serial number recorded with the data to assist in quality assurance. 
Additionally, the devices and corresponding phones were labeled 1, 2, and 3 (AirBeam 


001896105818, 001896105926, and 0018961061CE respectively) for easy identification in the 


field. 
x POWER BLUETOOTH LANYARD 
(A) BUTTON CONNECTION ATTACHMENT 
INDICATOR POINTS 





BATTERY 
CHARGING 
INDICATOR 


MICRO-USB 


INTAKE PORT 


BELT CLIP 
ATTACHMENT 


T -F - 
POINT emperature Sound Level - dB 


AirBeam-F Phone Microphone 


EXHAUST 


ADD-ON 
SENSOR PINS 





Figure 12. A) AirBeam diagram (Heimbinder & Besser, 2014) and B) smart phone app example. 
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Figure 13. The inside of the Shinyei PPD60PV (Heimbinder, 2013). 


Since temperature, relative humidity, wind speed and direction (Hemond & Fechner, 
2014; Nowak et al., 2014; Tong et al., 2015; Tai, Mickley, & Jacob, 2010), precipitation, and 
change in pressure over time have been shown to account for 30 percent of the daily variability 
in PM25in Southeastern U.S. (Tai et al., 2010), a handheld Kestrel 4000 pocket weather meter 
(Nielsen-Kellerman, Nelson, PA) was used to collect the ambient temperature (accuracy: +1 C; 
resolution: 0.1 °C), relative humidity (accuracy: +3%; resolution: 0.1%), and wind speed 
(accuracy: +0.1 m/s; resolution: 0.1 m/s) aligned with wind direction at each site in ten-minute 
intervals. The wind direction was obtained using a compass by noting direction at maximum 
wind speed. Field tests were performed during periods of no precipitation to protect the 
AirBeam units. 

2.2.3 AirBeam Equivalency Tests — Following best practices (Lodge, 1988) and EPA 
recommendations (Williams et al., 2004) any device used for monitoring air quality should be 
suitably calibrated. Based on the standard reference instrument method, the AirBeams were 


compared with the state's TEOM continuous monitor located at the Columbus airport. Due to 


instrumentation differences (the TEOM has omnidirectional intake, while the AirBeam units 
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have unidirectional intakes), data resolution differences (TEOM has |-hour resolution and 
AirBeams have 1-minute resolution), and the limited lifespan of the AirBeams (2-hour battery 
life), it was determined intra-unit correction (through equivalency tests conducted prior to and 
during field testing) would be the best standardization method. 

The three AirBeams were assessed for their equivalency across a range of PM2.5 
concentrations from 0 to 177 ug-m° before the start of the field tests. These units were assessed 
for equivalency at low ranges by running the devices over the course of three days in an 
undisturbed room, in which the room’s HVAC vents were sealed and entry into the room was 
limited to conducting the test. The Austin Air HealthMate Plus® (an air purifier capable of 
removing particulates larger than 0.3 microns) was used to reduce particulate matter within the 
room after small levels of smoke were allowed into the room through a small opening in the 
window to compare response to stimuli over time. The three units were equal distance from the 
opening in the window and the air purifier. 

A second indoor test was conducted to compare the units’ response to high particulate 
concentrations. This test was performed by burning a 160 g carbon fiber vinyl ester specimen in 
a muffle furnace at 600 °C for 1 hour 47 minutes. The AirBeams were positioned in the fume 
hood 1 m above the furnace smoke stack. The fume hood was allowed to run prior to and during 
the burning of the carbon sample to ensure steady airflow. During this equivalency test, unit 3 
disconnected from the phone app and did not record 22 minutes of data. These 22 minutes were 
removed from the statistical analysis for all units. 

The first two tests were conducted in controlled indoor environments. The last 


equivalency test was conducted outdoors to assess response of the units to stimuli without the 


ability to control environmental factors. The AirBeams were set up outside equal distances and 
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downwind from an outdoor wood burning stove. Smoke was allowed to escape from the top of 
the stove for 5 minutes, and then the fire was squelched in the stove. The units collected data for 
5 minutes before and after the smoke, and the data were assessed to determine if the units 
responded to the stimuli at the same time. 

The Pearson product-moment correlation coefficient was computed between the three 
AirBeams and the state TEOM. Model II linear regression was used to assess the relationship 
among the three AirBeams for all equivalency tests combined. As with time-series data, 
autocorrelation was an issue in the results of all three equivalency tests. The appropriate lag was 
determined and autocorrelation corrected for each individual test (indoor test 1 lag = 21 min, 
indoor test 2 lag = 4 min, outdoor test 1 lag = 2 min) before using one-way analysis of variance 
(ANOVA) to compare the fecniis of the three units. 

2.2.4 AirBeam Temperature and Relative Humidity — The AirBeam’s relative humidity 
and temperature sensors are good indicators to ensure the devices are not overheating or 
oversaturated. AirBeams are programmed to shut off at 100 percent humidity (Heimbinder & 
Besser, 2014) and the Shinyei PPD60PV performs best at temperatures of 0 to 45°C (Shinyei 
Technology Co., Ltd). Jiao et al. (2016) and Mukherjee et al. (2017) did not address the 
performance of the temperature and relative humidity sensors housed in the AirBeams. The 
output from these sensors were compared to Kestrel 4000 data obtained during field testing and 
to hourly data from the weather station at the Columbus Metropolitan Airport (WBAN: 93842) 
[gathered from the National Centers for Environmental Information (NCEI, 2017b)]. 

2.2.5 Field Test — During February 2017 (winter, leaf-free period), particulate 


concentrations were measured in one-hour sampling sessions at 15 study sites (Figure 14), on 


rain free days, along major Columbus roads (Interstate 185, Highway 80, Highway 27, and 
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Highway 280), near busy shopping centers, and close to restaurants that produce smoke (e.g. 
Burger King). The roads that run along the study sites are considered principle arterial highways 
and experience average daily traffic ranging from 27,000 vehicles (Manchester Parkway) to 
66,000 vehicles (Highway 80). As access to power was limited, the devices were operated for 


one hour to avoid exceeding the two-hour battery lifespan. 
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shaped arrangements (yellow square with x) indicated across Columbus, Georgia. 


Sites were categorized as dense tree buffers (n = 5), small tree lines (n = 6), and U-shaped 
tree arrangements (n = 4). At eleven sites more than one location was tested so each tree 
arrangement type had 10 sample locations. Barriers, like fences, walls, and steep drop-offs in 
elevation, limited ability for more than one sample location at four sites. For each sample 
location two AirBeam units monitored particulate concentrations, one within the tree stand and 
one in the adjacent open area (Figure 10). In total, 60 sample locations (30 pairs of tree stand and 


Open areas) were monitored. 
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The study sites were chosen based on proximity to major roads and busy shopping 
centers. The closeness to smoke producing restaurants was not accounted for in the original 
experimental design, and sites were not picked with this feature in mind. Each site needed a tree 
buffer with neighboring open area and a higher density of conifers than deciduous trees (as 
testing occurred during winter months). Sample locations picked within sites were based on 
distance from source, with first locations close to the road and successive locations set farther 
back as spaced allowed. Detailed pictures of each study site with marked test locations can be 
found in Appendix B. It should be noted that sites were not tested randomly as wind direction 
needed to correspond to the direction of particulate matter source and access to some sites 
(Haverty’s, Lazyboy, Colony Bank, and the three churches) was limited to specific days and 
times. 

At each site AirBeams were set at the height of 1.7 m, the average adult height in the U.S. 
(Ogden, Fryar, Carroll, & Flegal, 2004). All three units were started at the same location close to 
the street, in the open area, and collected data for at least five minutes. While positioned facing 
the direction of the road, a compass was used to determine the unit orientation, and units 
remained facing this direction for the entirety of each testing session. The two units with the 
most similar peak and average data were identified and used for subsequent testing. The two 
selected units were placed within the tree buffer and adjacent open area at the same distance 
from the road. The units sampled for five to ten minutes at this location. Units within the tree 
buffer and adjacent open area were moved back farther (if space allowed) for an additional five 


to ten minutes. All units were moved back to the start location for the final five minutes. At 


each location in the open area, the Kestrel 4000 was held at approximately 1.5 m and 
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temperature was allowed to equilibrate before wind speed, relative humidity, and temperature 
data were recorded. 

Aakash et al. (2017) recommends frequent recalibrations in an environment similar to 
study conditions when using low-cost portable devices. Therefore, the five-minute start and end 
PM2.s5 data were used to correct for AirBeam discrepancies. The units’ median values for the 
start/end period were calculated. A median ratio (unit 1: unit 3 and unit 2: unit 3) was then 
applied to the PMz.s data recorded at each location. Unit 3 was used as the standard for 
comparison testing as unit 3 was used at 27 of the sampling locations for tree/open concentration 
monitoring. (This matches the equivalency test results as unit 3 was found to vary between unit 
1 and unit 2). For the three locations in which unit 3 was not deployed, a median ratio of unit 1: 
unit 2 was used to correct baseline differences. While a high correlation was found during 
equivalency tests among the units, this added corrective measure was used to ensure the data of 
importance, the tree and open PM2s concentration data, aligned before final analysis. At three 
sites (the first site for each tree arrangement type), the units were corrected based on median 
ratio data for all sites and equivalency tests because the median ratio correction method 
(described above) was not employed until after testing was completed at these sites. 

Columbus hourly weather from the weather station at the Columbus Metropolitan Airport 
(WBAN: 93842) airport for February 2017 were obtained from the National Centers for 
Environmental Information (NCEI, 2017b). Columbus hourly PM2.s concentrations for February 
2017 were certified and obtained from the Georgia Environmental Protection Department - Air 
Branch Division. These data were used to assess if a relationship existed at city level between 


weather conditions and particulates. The daily averages were calculated and correlations 


examined between wind direction, wind speed, temperature, relative humidity, precipitation, and 
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Columbus PM2.s5 concentrations. Fort Benning controlled burn (Fort Benning's Smoke and 
Sound Archive, 2017) and regional agricultural fire (NESDIS, 2017) dates were gathered to 
determine if regional smoke impacted PM2.s and study sample locations’ PM2.s concentrations. 

The Columbus hourly PM2.s concentration were also used to assess whether study sample 
locations’ PMz2.s levels were influenced by localized sources. Each location’s peak PMz2.5 values 
less the city PM2.5 data were calculated. The results ranged from 2.0 to 92.6 u g-m™°. Values 70 
to 499 ug-m°® were labelled as high, 30 to 69 ug-m® as medium, and below 29 yg-m” as low 
sources. As a part of their Village Green Project, the EPA developed this scale for short-term air 
sensors in order to understand personal exposure to nearby air pollutants (Keating et al., 2016). 
This method matched what was experienced at the sites as the three high PM source locations 
were near local restaurants producing smoke and the one medium location occurred near a stop 
light with idling vehicles. These PM source level results were considered a random factor in the 
analysis as this was controlled for in the experimental design. The mean PM2:s for each location 
(tree vs open) and by type (dense, small, U-shaped) was used in an analysis of covariance 
(ANCOVA) with temperature, relative humidity, the wind direction versus the AirBeam unit 
orientation, and wind speed as covariates. IBM SPSS Statistics (Version 25.0) computer 
software was used for the statistical calculations (IBM Corp., 2017). 
2.3 Results 

2.3.1 AirBeam Equivalency Tests — Two of the three AirBeams had good correlation 
(Unit 1: r= 0.78, p = 0.066; Unit 2: r = 0.85, p = 0.032; Unit 3 r = 0.85, p = 0.03) to the state 
TEOM when PM? concentrations ranged from 0 to 10 ug-m°. The among unit equivalency 


tests yielded a total of 774 minutes (601 min for indoor test one, 85 min for indoor test two, and 


88 minutes for the outdoor test) of PM2.s data for each AirBeam. All three units showed a 
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significant positive relationship to each other (Figure 15; Model II Regression Units 1 & 2: Fiz? 
= 67471, r2 = 0.989, p < 0.001; Units 1 & 3: Fi,773 = 37019, r? = 0.980, p < 0.001; Units 2 & 3: 
F1,773 = 53262, r? = 0.986, p < 0.001). All three units were not statistically different from each 
other when comparing the mean PM2.5 among all three equivalency tests (ANOVA F2.278 = 
0.440, p = 0.645). Overall, units 2 and 3 differed least as compared to unit 1 (Units 2 & 3 Tukey 


HSD p = 0.991; Units 1 & 2 Tukey HSD p = 0.738; Units 1 & 3 Tukey HSD p = 0.661; Table 6). 
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Figure 15. Airbeam PM2.s equivalency test results showing model II linear regression 
relationship between A) unit 1 v. unit 2, B) unit 1 v. unit 3, and C) unit 2 v. unit 3. 1:1 line 
denoted by solid line for reference. 
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Figure 16. Equivalency test results PM2.s mean and 95% CI: A) indoor test one, D) indoor test 
two, C) outdoor test, and D) all tests combined. 


To better understand the relationship among the units, comparisons were made among the 
three separate equivalency test results to determine how the units performed with various 
particulate concentrations and indoors versus outdoors (Table 6). The three tests varied in PM2.5 
concentration ranges: indoor test one ranged from 1 to 10 u g-m”™, indoor test two ranged from 30 


to 177 ug-m®, and the outdoor test ranged from 17 to 155 u g-m-. No statistically significant 


difference was found between the average PM2.5 estimates for the three units for the first indoor 
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test (ANOVA F2,81 = 2.322, p = 0.105), the second indoor test (ANOVA F2,62 = 2.562, p = 
0.086), and the outdoor test (ANOVA F2,131 = 0.927, p = 0.398). Units 1 and 3 varied less at 
lower concentrations (4 percent difference) versus higher concentrations (20 percent difference). 
The pairwise relationship results were similar between the first indoor test and the outdoor test. 
Units 1 and 3 have a 4 percent mean difference for the outdoor test (Figure 16). 


Table 6. Equivalency test results for AirBeams' PM2.5 mean, 95% CI, and Tukey HSD p-values 
indicating no significant difference between PMz2.s concentration means of the three units. 


Indoor Test 1 Indoor Test 2 Outdoor Test Combined 
(n = 28) (= 21) (n = 44) (n = 93) 
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2.3.2 AirBeam Temperature and Relative Humidity — For average relative humidity, all 









three AirBeam units, the Kestrel and the City data were not significantly different (ANOVA F4,89 
= 2.425, p = 0.054; Table 7). The mean relative humidity for unit 3 was significantly different 
(at an a = 0.1 threshold) because it was 15 percentage points lower than Kestrel and City data 
and 7 percentage points lower than units 1 and 2. For temperature, a statistically significance 
difference was found (ANOVA F439 = 2.425, p = 0.001). Units 1, 2, and 3 showed mean 
temperatures that were 4 C and 5 C higher than the Kestrel and City, respectively (Figure 17). 


Table 7. Relative humidity and temperature results: Kestrel, City, and AirBeams’ mean, 95% CI, 
and Tukey HSD p-values. 


Relative Humidity (%) Temperature (°C) 
Kestrel City Unit1l Unit2  Unit3 | Kestrel City Unit 1 Unit 2 Unit 3 
Mean 45 45 37 37 30 21 20 25 =O 25 
+95%CI | +8.6 +9.0 +7.0 Sei +7.0 +1.8 10) +2.8 aa | +2.6 
Kestrel 1,000: 07553 0-534 0.074 0.951 0.090 0.102 0.080 
City 0.631 0.609 0.098 0.013 0.016 0.012 


Unit 1 1.000 0.785 1.000 1.000 
Unit 2 0.821 1.000 
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Figure 17. Hourly mean and 95% CI of A) relative humidity and B) temperature for Kestrel, 
City, and the three AirBeam units. 


2.3.3 Field Test — The temperature obtained using the Kestrel 4000 during the field 
testing period ranged from 12.1 to 28.6 ‘C with mean of 21.2 C, the relative humidity varied 
from 19.4 to 87.1 percent with mean of 51.1 percent, and the wind direction was variable at 
speeds of 0 to 5.1 m/s with a mean of 1.6 m/s. The temperature, relative humidity, wind 
direction in relation to unit orientation, and wind speed were assessed using linear regression and 
found to account for little variation in PM2:s across all locations (17 = 0.110, F4,23 = 0.712, p = 
0.592). Therefore, the Kestrel obtained weather data were not included in the rest of the 
analysis. 

The PM2>5 concentrations differed across PM source level (F252 = 41.635, p < 0.001). 
The tree stand arrangements also differed significantly (ANCOVA F2,52 = 3.939, p = 0.026). No 
difference in PM2.5 was found in open areas versus trees (0.9 percent difference) across all sites 


(ANCOVA F152 = 0.003, p = 0.956). Dense tree buffers differed as compared to small tree lines 


(43 percent difference; Tukey HSD p = 0.001) and U-shaped arrangements (31 percent 
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difference, Tukey HSD p = 0.018). However, small tree lines and U-shaped arrangements 
showed no differences in mean PM2.5 (Tukey HSD p = 0.633). Dense tree buffers had higher 
PM>s concentrations in trees versus open areas with a mean difference of 1.6 ug-m®. In 
contrast, small tree lines (mean difference of 0.1 wg-m™) and U-shaped (mean difference of 1.7 
ug-m™) arrangements had higher PM2.5 concentrations in open areas versus trees (Figure 18). 
This interaction was not statistically significant (F252 = 0.716, p = 0.493). Sample location 
details used in the analysis, including PM2.s corrected data and weather data, are located in 


Appendix B, Table 7. 
PEs 
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Figure 18. Comparison of the open and tree PM2.s concentrations mean and 95% CI within dense 
tree, small tree line, and U-shaped tree stands. 


The Columbus city weather and PM2:s data for February 2017 were also analyzed to 
better understand the city’s particulates in relation to atmospheric conditions. During February 
2017, the city experienced 20 days with particulate concentrations less than 12 g-m°, 7.4 days 
with 12.1 — 35.4 ug-m™, 12 hours with 35.5 — 55.4 ug-m™, and 2 hours with 55.5 — 150.4 ug-m°?. 
All 24-hour averages were below 12 ug-m°, except for 5 days (February 5, 6, 11, 14, and 18), 


where concentrations were between 12.1 — 35.4 ug-m°®, with the highest being 27.3 u g-m™ on 


February 18". 








5] 


Winds from the south or east bring fine particulate concentrations to the city during Fort 
Benning controlled burns. Fort Benning conducted controlled burns on February Ist, 11th, 17th, 
24th, and 25" when winds were predicted to be from the north and west (Fort Benning's Smoke 
and Sound Archive, 2017). Overnight between February 17" and 18" the wind direction shifted, 
and winds from the southeast brought smoke into the Columbus area. This phenomenon was 
observed through an increase in PM2:5 concentrations and the highest concentrations all month. 
Rain occurred later on the day of the 18" and reduced airborne particulates. An increase in 
PM2.s concentrations on February 11 could also be linked to Fort Benning prescribed burns. 
The remaining three 24-hour periods (February 5, 6, and 14) were most likely due to agricultural 
burning in other parts of the state (NESDIS, 2017). No field test days overlapped agricultural 
burning dates, and three test days took place on the same day as controlled burns (February cas 
17", and 24"). The city hourly PM25 ranged from 0.7 to 11.9 ug-m™® with a mean of 5.7 pg-m”° 
during the hours when field testing occurred. Regional agriculture and Fort Benning fires were 
not found to contribute to field test particulate levels. 

When Fort Benning prescribed burn data were removed, statistically significant 
correlations were found between city hourly PM2.s concentrations and the following city weather 
variables: wind speed (r = -0.253, p < 0.001), relative humidity (r = -0.080, P = 0.043), 
temperature (r = 0.134, p = 0.001), and precipitation (r = -0.093, p = 0.018). When daily 
averages for the month of February 2017 were calculated, the city PMz.s correlated positively 
with temperature (r = 0.460, p = 0.024) and negatively with relative humidity (r = -0.433, p= 
0.035). No correlation was found with daily city PM2.s and wind speed (r = 0.166, p = 0.438). 


These findings suggest the daily changes in Columbus city PM2.5 concentrations could be 


attributed to daily temperature and relative humidity fluctuations (Figure 19). 
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Figure 19. Average hourly data variability in Columbus February 2017 PM2.5 as compared to A) 
relative humidity (RH) and B) temperature. 


2.4 Discussion 

2.4.1 AirBeam Performance - The stability, reproducibility, and reliability are concerns 
when using low-cost PM sensors (Rai et al., 2017). The three AirBeam units were assessed for 
these three factors during the equivalency evaluations and throughout the field study. Stability 
refers to the performance of the sensors remaining constant over a period of several months (Rai 
et al., 2017). The stability of the units was determined based on the change in unit median ratio 
between units over the course of the research. The median ratio for all tests for unit 1: unit 3 was 
1.17 and for unit 2: unit 3 was 0.93. At the beginning of the equivalency testing period (October 
2016), the median ratios were 1.01 (unit 1:3) and 0.92 (unit 2:3). At the start of the field testing 
median ratios were 1.03 and 0.91 (February 2017), and at the end of the field study they were 
1.22 and 0.85, respectively. The change in the median ratio may be due in part to units not being 
cleaned during the testing period. The makers of the particle sensor inside the AirBeam suggest 
a lifespan of two to three years for the sensor if it is cleaned properly (Shinyei Technology Co., 


Ltd.). Jiao et al. (2016) found a change in response with “days of use” for the AirBeams they 


tested over 168 days. Given the small change in median ratio from beginning to end, the 
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AirBeams were found to be stable for the short time period of use. Longer testing sessions could 
prove problematic if particles collect on sensors over time. 

Reproducibility is the difference in measurements found between similar devices (Rai et 
al., 2017). As PMz2.s concentrations increased, the differences among the three units also 
increased. While at lower levels all three units did not differ greatly, and units 1 and 3 
corresponded most closely. Units 1 and 2 correlated best at higher levels. Correlations among 
the units were the same in the field study as with equivalency tests. Particulate levels were low 
at many study sites with average concentrations of 12ug/m? for all locations. Units 1 and 3 were 
found to be closer in mean and peaks, and, therefore, were used in tandem at half the study 
locations. As seen in the equivalency tests, the three units had the same response to stimuli and 
were found to have high reproducibility. These findings match the two studies previously 
conducted in regards to the high correlation found between AirBeam units (Jiao et al., 2016; 
Mukherjee et al., 2017). 

The relative humidity and temperature sensors within each unit were used as an 
additional means to determine sensor reliability. The relative humidity and temperature sensors 
housed within the units were used as a gauge to ensure the units would not shut off, as AirBeams 
are programed to shut off at 100 percent relative humidity (Heimbinder and Besser 2014) and do 
not operate outside the temperature range of 0 to 45 C (Shinyei Technology Co., Ltd). The 
AirBeams’ internal temperature and relative humidity were determined to not interfere with unit 
performance in the present study. The maximum relative humidity reading of all three units was 
80 percent recorded on unit 2. The maximum temperature reading was 38 “C and the minimum 


reading was 12 °C, both measured on unit 1. The comparison between the AirBeams, Kestrel, 


and City weather data suggest that these sensors should not be used when collecting local 
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weather data. For example, all three units’ mean temperatures were 17 percent higher than 
Kestrel and 22 percent higher than City mean temperatures. The AirBeams’ black casing may 
account for the higher average temperatures seen in each unit as compared with the Kestrel and 
City data. Using the units during the summer in Columbus, Georgia, could prove problematic as 
temperatures average 33 °C in July and August (NCEI, 2017a) leading to a greater probability of 
exceeding the acceptable AirBeam internal temperature range. 

The utility of the AirBeam PM25s sensors for extended deployments and over multiple 
seasons with varying temperatures has yet to be established. One additional problem is the need 
to tether the device to a cell phone with a data plan to record data in the field. However, the 
AirBeams’ high correlations, short-term stability, and reliability make these units legitimate 
sensors for the field studies. These units have limitations in data accuracy and usable 
temperature ranges, but the practicality of these devices is in their portability and ease of use. 
Field assistants were trained in less than five minutes to use the units and the accompanying cell 
phone app. The only complication occurred when field assistants did not save data properly. 
These features were found to be valuable when conducting the field study tests at multiple 
locations. 

2.4.2 Fine Particulate Matter and Trees — The relationship between PMz:5 and trees was 
statistically insignificant for all tree versus open sample pairs. While not significant, isolating by 
tree design found results similar to other studies. In this study, small tree lines had no impact on 
particulate levels as compared to surrounding open areas. This finding parallels the results of 
Hagler et al. (2012) that roadside buffers less than 10 m did not hinder particle transport beyond 


the vegetative barrier. Tong et al. (2015), Whitlow (2013), Liu et al. (2015) and Chen et al. 


(2015) found PM2.s concentrations were higher in dense tree buffers versus neighboring open 
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areas. While this study found no significance between all open versus dense tree buffer sample 
sites, it did find higher PM2.5 concentrations in dense field sample locations as compared with 
other tree configurations. Overall the small difference between open and tree concentrations 
when using active monitors is consistent with the observations of Setala et al. (2013). 

For dense tree arrangements, the trees trap nearby PMz2:s keeping particulates from 
leaving these tree barriers after being intercepted by the trees. Three locations used in this study 
did not have this relationship dynamic. For the Havertys and Lazyboy sites, PM2.s was 
essentially the same in the trees (m = 6.6 wg-m"*) and the open (m = 6.9 ug-m*). These locations 
were tested in the same sampling session because these sites are located near the same high 
traffic shopping center. The notable differences versus the other seven locations was the lower 
than anticipated traffic volumes, 3°C higher temperature, and PM2.s that was 12 ug-m™°? lower 
than the other dense field locations. It is not possible to determine whether the higher open 
PMz2.s concentrations is a confounding variable of these three locations, that, when removed, 
would make the overall dense tree buffer configuration results significant. Four of the dense 
field sites had a second location at farther distance from PM source. At each of these sites, a 
decrease in particulate concentration was measured with distance from PM source in both open 
and tree locations. 

One dense tree location, the Manchester bike park, had the highest overall particulate 
matter concentrations of the entire study (tree m = 39.8 wg-m™ and open m = 34.1 ug-m*). This 
location is located near a traffic stop, and cars and trucks were idling at the stop light with wind 
direction coming from the direction of this source during the sampling session. An increase in 


PM2.s is expected with idling vehicles, especially diesel trucks (Girard, 2014; Reff et al., 2009). 


The city PM2.s was 8.5 ug-m° during the test. Two other testing sessions (not used in statistical 
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analysis due to partial loss of data) at this location, saw the same relationship of higher 
particulates in the trees versus open areas at lower PM2:5 levels (5.8 and 5.4 ug-m™>). The 
sessions with data loss were conducted at low traffic periods, so this may account for the 
difference in particulate levels among sessions. 

Another dense tree site located on Williams Road near the I-185 exit 12 on/off ramp had 
relatively elevated particulate levels (average of tree m = 17.0 ug-m* and open m = 15.3 wg-m?). 
This site (Figure 10A) offered the perfect dense field set-up with dense trees to the east of a 
cleared open area and across the street from two gas stations. The slightly elevated PMz.s in the 
area was thought to be due to proximity to these gas stations and idling traffic. The winds were 
calm for the majority of the testing session with the exception of the start. The lower winds 
could also lead to a build-up of pollution in the area (Tai et al., 2010). 

The small tree line arrangement overall saw no statistical difference in tree PM2:s 
concentrations versus open PM2.s concentrations. The smaller number of trees in these 
arrangements are not able to noticeably reduce the fine particulate matter in the air. One site (the 
CSU softball field site) experienced high PM source levels, and open PM2.5 concentrations were 
2.9 ug-m™ higher than tree concentrations. This may be a feature of this location. Tall pine trees 
populated the site. Wind direction aligned to bring smoke from a local Burger King to this site. 
Smoke from a meat cooking restaurant is higher in the air column, and the tree tops should 
intercept some of the particulate pollution. The parts of Columbus that have small tree line 
arrangements are shopping centers with meat cooking restaurants. As discussed in Chapter 1, 


often smaller ornamental trees frequent these areas as compared with the trees found at this site. 


The CSU softball site, with additional testing, could serve as a case study example of how tall 
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trees within small tree line designs can reduce fine particulate concentrations, which may 
influence the types of trees planted near restaurants that produce smoke. 

While no relationship was found between field measured weather parameters and 
particulates, the city PM2.5 and city temperature, relative humidity, and wind speed were found to 
have a significant correlation. This weather-particulate interaction at the city level could help 
explain some of the findings at sample locations. The direction of relationship found between 
city PM2.s concentrations and city temperature, relative humidity, and wind speed are consistent 
with Tai et al. (2010) findings for Southeastern U.S. and point to organic carbon (OC) and 
elemental carbon (EC) in the atmosphere. OC and EC are mainly caused by combustion of fossil 
fuels, which is consistent with Fort Benning controlled burns, vehicles, and smoke from 
restaurants as sources of particulates. Additionally, the daily decrease of particulates from 
morning to afternoon might explain the higher levels of particulate matter found at Cascade Hills 
Church (small tree line site) during the only sampling session that took place in the morning. All 
three locations sampled at the church averaged 12 ug-m™ the morning tested with relatively high 
traffic conditions (83 vehicles per minute), but two previous tests conducted in the afternoon 
measured particulate levels below 3 ug-m™ with higher traffic conditions (99 vehicles per 
minute). 

The U-shaped tree arrangement also had no statistically significant difference between 
tree and open area particulate concentrations. At all PM source levels open area concentrations 
were higher by 1.7 ug-m™* when compared with tree concentrations. Not all study sites 
constituted ideal U-shaped tree stand arrangements, making city-wide generalizations difficult. 


However, the All Saints Presbyterian Church (Figure 10C) could be considered an almost perfect 


U-shaped tree stand with idling cars and diesel trucks at the opening of the U being the main 
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source of particle pollution. The average PM2:s at this site was low, but high winds (the highest 
recorded throughout the entire field study at 5.1 m/s) from the direction of the road brought an 
increase in particulate levels to the open area as compared to the trees. While sampling at the 
second location, the winds increased from 2.9 to 11 m/s. Unit 2 was left at the start/end location 
approximately 35 m from the road. Unit 1 was positioned 85 meters from the road in the open 
area, and unit 3 was the same distance in the tree line. With the increase in winds the particulate 
level also increased. This was seen with particulate levels peaking in series four times, first at 
unit 2 and 25 to 37 seconds later at unit 1. This dynamic of flowing through the opening of the U 
and not the trees highlights the impact the right wind direction and wind speed can have on 
particulate levels in this tree arrangement. 

Two other U-shaped sites (Colony Bank and the corner of University Avenue and 
Manchester Expressway) were located across the street from restaurants that produced smoke 
during testing sessions. Wind direction was from the direction of these sources, and, 
consequently, these sites experienced high PM2.s concentrations. The Colony Bank site is also 
located near a shopping center parking lot. The U-shaped opening points towards this parking 
lot, while a small tree line exists directly opposite the restaurant. The Colony Bank site had the 
same average PM2s levels (14.6 ug-m’) in the trees and open areas with a slightly elevated level 
in the trees (22.1 ug-m™) as compared to the open (21.9 ug-m~) when smoke was present. 
Conversely, at the corner of University Avenue and Manchester Expressway the U opens 
towards the restaurant. The average particulate concentrations were higher in the trees (12.9 ug- 
m™) versus the open (8.1 ug-m™) being 8.5 ug-m> higher in the trees over the open area when 


smoke was present. A bike path runs through the U-shaped opening at this site, meaning higher 


particulate levels are experienced by people using this path for recreation when the restaurant is 
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cooking meat. The difference between the two sites demonstrates how the right (or wrong) 
alignment of trees to PM source impacts particulate levels. The city of Columbus has several of 
these U-shaped tree designs because often only trees necessary for development are cleared in 
order to save tree canopy. This tree arrangement becomes problematic when it is located in areas 
where people frequent, like parks, and a PM source is near. 

Every site, even those closely located (i.e. Manchester bike park and the corner of 
University Avenue and Manchester Expressway), has different localized PM sources that 
contribute to particulate levels. As discussed in methods and results sections, these localized 
sources must be controlled for in order to compare tree stands across the city. It is important to 
note these PM sources beyond controlling for background though. A study conducted in 
southeastern United States cities found wood combustion made up 25 to 66 percent, diesel 
exhaust 14 to 30 percent, meat cooking operations 5 to 12 percent, and vehicle exhaust 0 to 10 
percent of the OC PMz2:s concentrations (Zheng et al., 2002). The portion of fine particulate 
matter caused by vehicles in Atlanta, Georgia, has decreased due to vehicular emission 
regulations (Vijayaraghavan et al., 2012). The low wind speeds in Columbus cause vehicle 
particulate pollution to remain in the vicinity of the roads. This concept could be seen with low 
particulate levels at most locations near areas of high traffic (sample locations were greater than 
15 m from the road), except those mentioned already as being located near smoke producing 
restaurants or heavy idling traffic with diesel trucks and cars. 

Smoke producing restaurants were the sources of high PM in this study. Smoke from 
restaurants were found to pass through sampling areas within minutes. These higher 


concentrations plumes (73 to 93 ug-m™) can impact people in sensitive groups such as those with 


respiratory issues (EPA, 2016b). A study conducted in Los Angeles found meat cooking 

















60 


establishments contribute 21 percent of the OC PM2:s concentrations in the city (Rogge, 
Hildemann, Mazurek, Cass, & Simoneit, 1991). Reducing pollution at the source is the best way 
to combat it, but cities have shown little will to regulate restaurant emissions (Murphy, 2015; 
Chaudhury, 2015). Taller trees can assist in blocking the spread of smoke from these point 
sources, as seen at the Colony Bank and CSU softball field sites. Future studies should focus on 
this dynamic looking at height of trees near smoke producing restaurants as well as distance to 
source. 

This field study test had limitations. The testing took place for one month, only 
encompassing one season of the year. The study, while city wide, was on a small scale based on 
the number of locations visited. The small sample size limits the ability to apply findings 
beyond specific sites tested. The use of the Kestrel 4000 and its limitations may be the main 
reason for the insignificant statistical relationships seen between PM2.s and weather conditions 
measured at each site. The Kestrel 4000 is not capable of determining wind direction. Wind 
speed and direction were variable and Kestrel sampling was not continuous, rather a sample 
point method was employed. Continuous weather monitoring with similar sample resolution as 
the AirBeam is needed to assess if localized weather influenced PM2.5 concentrations. 

The time of year offers some complications as deciduous trees had shed their leaves 
before this study took place. Leaf absorption of ultrafine particulates is limited (Hemond & 
Fechner, 2014). Fine and ultrafine particulates settle on leaves through deposition. PM2.5 levels 
on leaves are lower than PMio due to gravitation deposition properties (Beckett et al., 2000; 
Freer-Smith, 2005; Seb¢ et al., 2012). Conifers are better at capturing particulate matter due to 


leaves having a waxy coating, (Szb¢ et al., 2012), high leaf area index, and no annual loss (Yang 


et al., 2015). Broadleaf trees have the second greatest capacity to capture airborne particles 
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(Beckett et al., 2000; Yang et al., 2015). All study sites had conifers trees (site pictures in 
Appendix B show sites during leaf-off season). Focusing on locations that had more conifer tree 
species during the winter makes findings specific to study locations during the one season. A 
Beijing study comparing forest PM2.5 concentrations to open area concentrations found the same 
higher concentrations in forests during leaf-off periods (Liu et al., 2015), while Cai et al. (2017) 
found higher deposition levels in urban settings in winter months. The composition and main 
sources of fine particulate matter can change throughout the year, with more wood combustion in 
the winter and higher biogenic VOC in the summer (Tai et al., 2010; Malm, Schichtel, Pitchford, 
Ashbaugh, & Eldred, 2004). These changes could have an impact on tree-particulate interactions 
throughout the year in addition to leaf-on vs leaf-off differences (Cai et al., 2017). 

Additional field study tests need to be conducted to determine the appropriate level of 
tree services in reducing PM2:s5 taking into account various localized particulate sources. The 
insignificant statistical findings between open and tree PM2.5 concentrations point to the low 
ability of conifer trees at these study locations to trap particulates during the winter season. A 
larger sample size, across multiple seasons will help in determining if similar findings are 
significant annually and city-wide for the city of Columbus. Dense tree barriers may reduce 
PM2.s concentrations in other seasons. Additionally, taller trees may assist in the reduction of 
airborne smoke particulates from nearby restaurants. The results of this study highlight the need 
to focus on various tree configurations in relation and distance to different particulate sources 
when considering utilizing trees as a deterrent to saeteulare pollution. Low-cost, portable 
sensors, like the AirBeam, can aide in determining neighborhoods with higher relative PM2.5 


concentrations and identify sources, as well as, assist in determining appropriate tree design and 


placement to reduce pollution. 
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DISCUSSION 

When looking at air quality benefits, the PM2.5 results found using the i-Tree model in 
Chapter | should be discussed with respect to the results of the PM2:s field tests conducted in 
Chapter 2. The field tests (i.e. Chapter 2) can better characterize the interactions between trees 
and PMz.s5 on a site by site level and facilitate generalizations regarding if the PM removal 
associated monetary savings are valid for the Columbus area. The field study documented 
higher particulate levels in treed versus adjacent open areas in seven of the ten dense field 
sample locations tested. On average, the observed difference was not great (1.6 u g-m~). The 
greatest variation in PM2.s was measured at the Manchester bike park, the area with highest 
average particulate levels for the whole study. 

Most of the northern portion of Columbus consists of dense fields of tree. Census tracts 
103.02 (Bradley Park area, south of 102.01), 33.01 (area northwest of the I-185, highway 280 
intersection), 105.02, and 105.01 (to the west of 108.02) are less developed with higher tree 
canopy percentages and dense tree buffers are the main type of tree arrangement (See Figure 8 or 
Appendix A, Figure 2 to locate tracts). Assuming the i-Tree Tool removal rates are accurate, 
these nine tracts contain 14,516 ha (35,871 acres) of canopy that could potentially remove 15 
tons of particulates (or 18 percent of the original overall city removal of this pollutant) during 
Columbus winters (i-Tree Tool removal rate was adjusted to account for the field study being 
conducted during leaf-off period). The notion this removal rate is valid for these areas hinges 
not only on dense tree stands trapping particulates, but also on particulate pollution being the 
same in these nine tracts as it is for the whole city. This rate depends of pollution levels gathered 
at the Columbus Airport. Witlow (2009) argues localized particulate concentrations differ from 


those detected by regional monitors. These nine tract areas are less populated with less traffic 
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and meat cooking restaurants, and, therefore, the fine particulate levels may not be as high as 
those near the airport. Regional fires may contribute to particulate pollution in these mainly 
northern tracts with the right wind direction. However, based on the tested conducted during 
controlled burns at increasing distance from Fort Benning by Baumann (2005) and Liu (2010), 
the smoke from Fort Benning controlled burns will likely disperse before significantly impacting 
these areas. 

The remaining 44 census tracts contain residential, business, and shopping neighborhoods 
with various tree arrangements. PM2.s removal by trees is not as accurate in these areas using the 
i-Tree Tool. Additionally, if the results from the field study hold, U-shaped tree stands have 
open areas with elevated particulate pollution, and small tree stands would not impact this 
pollution. Therefore, the areas of the city with the greatest population would not observe trees 
reducing particulate levels. 

In addition to limitations previously discussed, the applicability of the field test results 
relative to the i-Tree PM2.5 removal rates is questionable. The field study took place over the 
course of one month (February), while the i-Tree removal rates are annual rates. Simply 
adjusting the rate to cover one season, as above, does not account for the change in particulate 
levels, sources, and trees across all seasons. While Columbus, GA, experiences similar wind 
conditions throughout the year, other seasonal factors such as tree leaf off, weather, and fine 
particulate matter variations were not taken into account in this research. One big difference 
between summer and winter seasons is the existence of more leaves on trees to intercept 
particulates. Trees emit more volatile organic carbons during the summer, which increases 


ozone and can lead to eventual increase in particulates within and around trees (Yang et al., 
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2015). These seasonal dynamics could change the interactions observed between trees and 
particulates in the field study. 

Neither studies’ results unequivocally conclude that trees abate PM2.s in Columbus. 
Additional research is needed to assess the effectiveness of trees for reducing particulates 
specifically as tree planting activities relate to particulate reduction. Researchers have argued 
that planting trees solely for the purposes of improved health from reduced particulates is in 
“vain”, and that government funds should be used to reduce pollutants at their source (Whitlow 
et al., 2014). While this argument is valid, as discussed in Chapter 2, not all pollutant sources 
are regulated at the source. A telling example is the lack of desire to control pollutants from 
meat cooking restaurants (Murphy, 2015). Trees also help cities in other ways, like cooling air 
temperatures, reducing storm water runoff, and improving health not related to air quality (Pataki 
et al., 2011). Therefore, the aim of tree planting should be to provide maximum total benefits of 
the services offered by trees as a whole and not simply their ability to remove particulate 
pollution. 

In this research, the ability of trees to reduce air pollutants was assessed using high 
spectral analysis that quantified tree canopy and its associated benefits. The overall the canopy 
for Columbus, Georgia, at 52 percent, meets the criterion set by the American Forests Urban 
Forest Program for ideal urban canopy cover (Leahy, 2017), but the variations in percent cover 
across the city leaves the impervious downtown, business, and shopping center areas lacking in 
good canopy cover. Urban canopy cover recommendations are made so cities can benefit from 
the ecosystem services trees provide, but simply adding canopy does not mean these benefits are 
fully utilized. Tree placement, tree type, and tree design need to be considered, and the latter 


often is not when considering urban vegetation plans. 
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The city’s high tree canopy cover is estimated to remove 1,900 tons of criterion air 
pollutants and sequesters 282,000 tons of carbon dioxide annually. The high spectral imagery 
analysis highlighted large tree canopy and air quality benefit disparities over time across the city 
of Columbus. Areas of highest removal of gaseous air pollutants is dependent on location of 
trees. These are the northern sections of the city, which also have fewer air pollutant sources. 
Higher pollution and a lower number of trees in more urban areas of the city (downtown and 
shopping centers) lead to lower pollution removal. 

This research utilized a low-cost, portable particulate sensor to analyze the interactions 
between fine particulate matter and tree stand designs. AirBeams are affordable ($250/unit) and 
easy to use. Accurate, more expensive equipment, is not feasible for studies of this scale and 
length or reasonable for citizen use. The three units tested in this study effectively measured 
PMz2.s variations at multiple sample locations. The unestablished stability of the AirBeam over 
extended periods and its temperature restraints limits usability for longer testing periods. 
AirBeams and other portable sensors allow simple, city-wide field studies to be performed. Use 
of more affordable sensors leads to more measurements by more people, which in turn yields big 
data with incredible potential. Open access to big data allows for new possibilities in 
understanding the environment. The possibilities, given advancements in portable sensors, are 
wide, and can be very valuable in understanding air quality as it relates to many aspects of an 
urban setting at localized levels. 

PM2.s has more complex interactions with trees than other air pollutants and removal is 
dependent on local PM sources, weather conditions, and tree design. Small tree lines have no 
discernable impact on PM2.5 concentrations and dense tree buffers trap PM2.5 resulting in slightly 


higher tree particulate concentrations as compared to open areas. U-shaped tree stand 
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interactions with particulates depended on location of the open area within the tree stand in 
relation to notable PM sources. Overall, in the winter, trees had little impact on particulate 
concentrations as compared with open areas. While the city of Columbus has, on average, low 
wind speeds, wind direction played a key role in particulates reaching sampling locations. 
Future tree plantings and removal should take note of areas with lower tree canopy as well as 
paying attention to tree arrangement and proximity to PM sources to better assist in the removal 
of air pollutants. Also, as the dense tree buffer arrangements trap pollution particles, the clearing 
of fields of trees should be seen as impeding the removal of PM2:s5 along with other air pollutants. 
Given the limitations of the study conducted, future research is needed to better 
understand the relationship between tree stand arrangements and fine particulate matter 
involving more sample locations across multiple seasons. Research should also focus on 
alternative fine particulate sources in addition to that from vehicles, like restaurants that produce 
smoke. Using portable monitoring devices to assess smoke fallout and interception by placing 
sensors in trees of varying height would be useful in determining tree height effect on local 
pollution produced by restaurants. Research is lacking in this area, and more portable sensors 


allows for more methods to asses these interactions. 
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Table 3. 2015 iso cluster values and reclassified values for 4 clipped sections. 
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Legend 
Difference 2005 to 2010 


RR 18 to 12 
a 11 to-s 


2 * Miles A 


Esri, HERE, DeLorme, Nepmylindis, © OpenSteetiiap contributors, and the GIS 
us 2 community 


Legend 
Difference 2010 to 2015 


a 6 to -2 


: * Miles A 


Esri, HERE, DeLorme, Msapmylndis, © CpenStreetMap contributors, and the GIS 
us er community 


Figure 1. City of Columbus tree canopy change by census tract between: A) 2005 and 2010 and 
B) 2010 and 2015. Dark green represents canopy gains and light green represents loss in canopy 
Over time. 
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Figure 2. Muscogee County, Georgia, census tract reference map — 2010 census (from 
census. gov). 
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Table 5. City of Columbus census tract percent tree canop 
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Table 6. City of Columbus census tract air quality benefits. 


Air Air 
NAME Area | Tree Pollutant CO2seq Pollutant 
(ha) | (ha) Removal (kg/yr) Removal 

(kg/yr) (kg/ha) 
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APPENDIX B — PM2.5 FIELD STUDY SUPPLEMENTAL DATA 
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Legend 

PM_ Median 
A 1.4- 10.0 
@ 10.1-20.0 


@ 201-308 


ces : Etri,.Delorme, USGS, NPS. Esri, HERE, DeLorme, Mapmylndis, 
Ss 


nSteetiiap contributors. ahd the GIS us er community 


Legend 


PM_ Median 
A 1.4- 10.0 
@ 10.1 - 20.0 


@ 201-398 


Sources: Esri, DeLorme, USGS, NPS, Esri, HERE, DeLorme, Mapmyindia, ©, 
OpenStreetMsp contributors, snd the GIS us er community 





Figure 3. A) Variation in average particulate concentrations by sample location. B) Zoomed view 
to highlight locations closely clustered. 








105 


Field Sites AirBeam Tests Metadata 

The following figures and metadata provide details on each study site. Figures highlight test 
locations, tree and open field, with white markers (numbers indicate which AirBeam unit was 
used at that place, i.e. #1, 2 or 3). The start and end location, equivalency check, of all AirBeams 


is represented with a yellow marker. 


46 Street River Walk Access Parking Lot 2-1-17 


When: 2-1-17 16:00-17:15 EDT 

Who: K. Youngquist and Care Bacon 

Where: 4th Street River walk access parking lot 

What: 4th street has a thin Ime of trees near highway 280. Test in tree ine compared to open 
grass near parking lot. Units 1 and 3 used based on equivalency tests. 

How: Unit 1 placed in treeline and unit 3 in field next to parking lot. 

Notes: 26 trees bigger than 3 inch diameter and less than 5 inch diamter in lne of sight. Road 


is higher elevation than trees and field. 


Hypothesis: Thin tree line will not create enough of a fence to reduce particulate matter farther 
from the road as compared with open parking lot without trees. 


Particulate Notes 
Time Note 
16:25-16:37 2 cars idle in parking lot 


e Smell of smoke. Fort Benning perscribed burn earlier (around noon). Wind shifted from 
16:54-17:01 out of SW 
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Weather 


Weather.com: 






Time 
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Airbeam Location Data: 
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Figure 5. The River Walk parking lot has a thin tree buffer 





djacent to highway 280. 
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Columbus State University ROTC 2-217 


When: 2-2-17 15:30-17-52 EDT 

Who: K. Youngquist and Care Bacon 

Where: CSU Lindsey Creek Road/ROTC Field 

What: CSU has dense line of trees near J-185. Test in tree ine compared to open area near sign. Units 1 
and 3 used based on equivalency tests. 

How: Unit 1 placed in treeline and unit 3 in open. 

Notes: Road is 10 feet lower elevation than trees and field. 

Hypothesis: In the winter, dense tree line will have higher particulate matter level as compared with open area 
without trees. 

Particulate Notes 


Campus police smoking near bridge. Smelled worse in trees. 
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Ambient Weather Data (Kestrel 4000): 
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Figure 5. Columbus State College has a big field of trees adjacent to an open space. 








109 


Manchester Expressway Park and Ride Bike Park 2-3-17 


When: 2/3/2017 16:30-17:18 EDT 
Who: K. Youngquist and Trevor Gundberg 
Where: Manchester Expressway Park and Ride Bike Park, 3690 Manchester Expy. Columbus, GA 31909 
What: The bike park has a dense patch of trees surrounding the parking lot and playground on all 
sides except the north entrance to parking lot. Test beyond tree line compared to parking lot. 
How: Unit #1 in open and #3 in tree line. Moved two airbeams at similar distances from road. behind 
trees and the other in parking lot. 
Notes: Broke pencils and pens, cut session short. Only one other car in parking lot due to cold, windy weather. 


Hypothesis: In the winter, the dense tree line will create a fence, reducing particulate matter farther from 
the road as compared with open parking lot without trees. 


Particulate Notes 
Area PM2?.5 via GA EPD Ag Branch: 





Weather 
Weather.com: Cloudy the whole time 
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AFF. 


ster Bike Park 2-3-17 


—_ 
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Figure 6. The Bike Park has U-shaped tree canopy with an open area in the center adjacent to 


Manchester Expressway. 


All Saints Presbyterian Church 2-9-17 


When: 2-9-17 16:30-17:30 

Who: K. Youngquist (Alone) 

Where: All Saints Presbyterian Church, 7170 Beaver Run Rd, Midland, GA 31820 

What: All Saints has a dense patch of trees surrounding the parking lot. Test in trees and 
beyond tree line compared to parking lot. 

How: Start/end three airbeams at distance from road. Pick two with closest averages and 
peaks. Leave third at start location. Move two comparable airbeams at similar 
distances from road, one in trees and the other in parking lot. 

Hypothesis: In the winter, trees will create a fence, reducing particulate matter farther from the 
road in trees as compared with open parking lot without trees. 

Particulate Notes 

Tine 


16:39 | Truck idling and all airbeams particulate count increased 


Car Data: 
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Weather 
Weather.com: 
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Figure 7. All Saints Presbyterian Church location represents a U-shaped tree arrangement with 
an open field in the center adjacent to highway 80. 
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Cascade Hills Church 2-10-17 


When: 2/10/2017 16:20-17:00 EDT 

Who: K. Younggquist and Trevor Gundberg 

Where: Cascade Hills Church, 34th Street, Columbus, GA 31904 

What: Cascade has a thin hne of trees east of the church building. Test beyond thin tree line compared to 
parking lot. 

How: Start'end three airbeams at fence boarder facing highway. Pick two with closest averages and peaks. 
Leave third at start location. Move two comparable airbeams at same distances from road, behind 
trees and the other in parking lot. 

Notes: Started second location, but church event caused early end. End equivalency test not performed due 


to manager informing us it was time to leave. 
Hypothesis: In the winter. trees will create a fence, reducing particulate matter farther from the road as compared 
with open parking lot without trees. 
Particulate Notes 
Car Data: 
2 Cars (Trevor | Cars (Kristin 
Time 
Count) Count) 
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Ambient Weather Data (Kestrel 4000): 


Wind Wind Wind Speed .«, | Dew Point | WetBulb | Pressure 
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Airbeam Location Data (My iphone): 
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2 5 a 
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Figure 8. To the east of Cascade Hills Church, the tr tree buffer is small aCe to MEER 80. 


Manchester Expressway Park and Ride Bike Park 2-13-17 


When: 2/13/2017 17:00-18:00 ET 

Who: K_ Youngquist and Trevor Gundberg 

Where: Manchester Expressway Park and Ride Bike Park 3690 Manchester Expy. Columbus, GA 31909 

What: The bike park has a dense patch of trees surrounding the parking lot and playground on all sides except the north 
entrance to parking lot. Test beyond tree ine compared to parking lot. 

How: Start'end three airbeams at distance 70 ft from road in grass north of parking lot. Pick two with closest averages 
and peaks. Leave third at start location. Move two comparable airbeams at similar distances from road, behind 
trees and the other in parking lot. 

Notes: 1 and 2 were not set to record until 5:21pm and 5:34pm respectively. Closest airbeams based on averages and 


peaks during start will be used moving forward (as was done on west side of parking lot) and not previous 
equivalency tests. 


Hypothesis: In the winter. the dense tree line will create a fence, reducing particulate matter farther from the road along tree line 
as compared with open parking lot without trees. 


Particulate Notes 


Car Data: 
Cars (Trevor | Cars (Kristin hieames Cars/Min 
Count) Count) 
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Area PM2.5 via GA EPD At Branch: 





Weather 
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Figure 9. The Bike Park has U-shape tree canopy with an open area in the center adjacent to 
Manchester Expressway. 
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Cascade Hills Church 2-15-17 


When: 2/15/2017 5-5:10pm 

Who: K. Youngquist and Trevor Gundberg 

Where: Cascade Hills Church, 54th Street, Columbus, GA 31904 

What: East of the church building, Cascade has a newly cleared openly in tree line. 

How: Started three airbeams at fence boarder facing highway. Moved all back same spot. 
Notes: Only did equivalency tests as church members started arriving. 


Weather.com: 
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Figure 10. To ho east of Cascade Hills Church, the tr hee parter is small Tent to RishWey 80. 
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Cunningham Center 2-15-17 


When: 2/15/2017 17:30-18:15 EDT 

Who: K. Youngquist and Trevor Gundberg 

Where: Cunningham Center, CSU, 3100 Gentian Blvd, Columbus, GA 31907 

What: The Cunningham Center has a thin patch of trees lining the street and part of the 


parking lot. Test tree line compared to parking lot. 


How: Start'end atrbeams at distance 30f from road in grass north of parking lot near Cunningham sign. 
Pick two with closest averages and peaks. Leave third at start location. Move two comparable 
airbeams at same distances from road_ behind trees and the other in parking lot. 


Notes: 3 was not set to record until 5:36pm. 3 fell over at 5:56 while being moved. 
Hypothesis: In the winter, the small tree line will not impact particulate matter as compared with 
open parking lot without trees. 


Particulate Notes: 
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Google Earth | ae) OA Se : 
Figure 11. Cunningham Center has thin tree line adjacent to open parking lot. 
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Manchester Expressway Park and Ride Bike Park 2-17-17 


When: 2/17/2017 16:00-18:00 EDT 

Who: K. Youngquist and Care Bacon 

Where: Manchester Expressway Park and Ride Bike Park, 3690 Manchester Expy, Columbus, GA 31909 
What: The bike park has a dense patch of trees surrounding the parking lot and playground on all sides except 


the north entrance to parking lot. Test beyond tree line compared to parking lot. 


How: Start'end three airbeams close to road in grass at west corner of park. Pick two with 
closest averages and peaks. Leave third at start location. Move two comparable airbeams 
at similar distances from road, within trees and along path. 


Notes: Testing started at 4:21pm, but airbeam 3 data (in the tree line) was not saved. So analysis 
can not be conducted. 

Hypothesis: § The trees will create a fence, increasing particulate matter in the tree line as compared to 
open parking lot. 

Particulate Notes 


Area PM?.5 wia GA EPD At Branch: 
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Google Earth ees ot 
Figure 12. The Bike Park has U-shape tree canopy with an open area in the center adjacent to 





Manchester Expressway. 


Corner of University and Manchester 


When: 2/17/2017 17:00-18:00 EDT 

Who: K. Youngquist and Care Bacon 

Where: Fall Line Bike Path - Comer of University and Manchester 

What: The bike path has a tunnel of trees surrounding the path. Test on path and in tree line at 
distance from street. 

How: Start'end three airbeams 10ft from road in grass at comer. Pick two with closest averages 
and peaks. Leave third at start location. Mowe two comparable airbeams at similar 
distances from road, within trees and along path in open area. 

Notes: Smoke from Burger King across the street started around 17:34 and was picked up by units 2 and 3 
farther from sources as compared with unit 1 closest to source at start/end location. 

Hypothesis: The trees will create a fence, reducing particulate matter in the tree line while increasing 
it within the tunnel created by the trees as move away from the road. 

Particulate Notes 

Car Data: 
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Figure 13. Trees create U-shape/tunnel arrangement around bike path at corner of University and 


Manchester. 











Colony Bank 2-19-17 


When: 2/19/2017 12:30-13:30 EDT 

Who: K. Youngquist 

Where: Colony Bank, 1531 Bradley Park Dr, Columbus, GA 31904 

What: The bank has U- shape arrangement of trees lining street at Bradley Park Drive and at the 


back of the bank between the parking lot and the highway 80 on ramp. Test at small opening 
in tree ine and in tree line at distance from street. 

How: Start/end three airbeams 15ft from road in grass at SE corner of bank lot. Pick two with 
closest averages and peaks. Leave third at start location. Move two comparable airbeams at 
similar distances from road, within trees and along opening in tree line. 

Notes: Winds from NNW not in right direction for most traffic idling at stoplights. but only had 
permission for two days. So moved forward. 

Hypothesis: Wind direction will have greater impact reducing tree baracade effect. The dense tree line 
will not reduce particulate matter farther from the road more as compared with open parking 
due to wind direction. 


Particulate Notes 
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Airbeam Locations: 
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Figure 14. View of Colony Bank and surrounding shopping center area. 
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Figure 15. Trees create U-shape tree arrangement around Colony Bank parking lot. 
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Haverty’s and Lazyboy 2-19-17 


When: 2/19/2017 14:57-15:55 EDT 

Who: K. Younggquist and Will Kiourtsis 

Where: Havertys and Lazeboy, 5555 Whittlesey Blvd #1000, Columbus, GA 31909 

What: Line of trees at fence separating back of Havertys store from exit/on ramp to highway 


80 at Veterans Parkway. Traffic sits at street ight waiting to turn onto Veterans. 
Possible PM build-up at ight. Similar situation at street near Lazyboy. 

How: Start/end three airbeams at opening in tree line along fence. Pick two with closest 
averages and peaks. Leave third at start location. Move two comparable airbeams one 
in trees and two in openings. 

Hypothesis: The tree dense line will trap particulate matter increasing levels in trees as compared 
to open areas. 


Particulate Notes 
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Figure 17. Back side of Lazyboy has dense field of trees next to grass opening. 
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Colony Bank 2-20-19 


When: 2/20/2017 17:15-18:15 EDT 

Who: K. Youngquist and Trevor Gundberg 

Where: Colony Bank, 1581 Bradley Park Dr, Columbus, GA 31904 

What: The bank has U-shape arrangement of trees lining street at Bradley Park Drive and at the 


back of the bank between the parking lot and the highway 80 on ramp. Test at small 
opening in tree ine and in tree ine at distance from street. 

How: Start'end three airbeams 15ft from road in grass at SE corner of bank lot. Pick two with 
closest averages and peaks. Leave third at start location. Move two comparable airbeams 
at similar distances from road, within trees and along opening in tree line. 

Notes: Smoke from Burger King during testing. 

Hypothesis: The dense tree Ime will create a fence, reducing particulate matter farther from the road 
more as compared with open parking lot without trees. 
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Airbeam Location: 
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Figure 18. Trees create U-shape tree arrangement around Colony Bank parking lot. 


Cascade Hills Church 2-23-17 


When: 2/23/2017 7:35-9:05 EDT 

Who: K. Youngquist and Dalton Peters 

Where: Cascade Hills Church. 54th Street, Columbus, GA 31904 

What: Cascade has a line of trees near the entrance to parking lot and a second line 
of trees past the church building. Test in tree line compared to parking lot. 

How: Start'end three airbeams at fence boarder facing highway. Pick two with 


closest averages and peaks. Leave third at start location. Move two 
comparable airbeams at similar distances from road, behind trees and the other 
in parking lot. 

Hypothesis: The highway traffic is not close enough to impact PM levels at the church. 
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Area PME 5 via GA EPD Aw Branch: 
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Figure 19. Cascade Hills Church the tree buffer runs parallel to highway 80. 
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Columbus State University Softball Field 2-23-17 


When: 2/23/2017 16:11-17:15 EDT 


Who: K. Youngquist and Kiara Mills 
Where: North of CSU Softball Field, 3100 Gentian Blvd, Columbus, GA 31907 
What: To the north of the CSU softball field, pine trees line the street. Test within tree 


line compared to open field. 

How: Start/end three airbeams in grass north of parking lot near end of north 
Cunningham building. Pick two with closest averages and peaks. Leave third at 
start location. Move two comparable airbeams at similar distances from road 
within trees and the other in parking lot. 

Notes: Airbeam 2 had trouble connecting. So used 1 and 3 for test. Ended longer to 
record Airbeam 2 at higher PM levels. 

Hypothesis: The tree tops are not full enough due to pruning to reduce particulate matter 
more as compared with open area without trees. 
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Airbeam Location: 
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Figure 20. North of CSU softball field is a small patch of trees. 


St. Mary’s United Methodist Church 2-24-17 





When: 2/24/2017 1552-17-06 EDT 

Who: K. Youngquist and Kiara Mills 

Where: St. Mary’s Road UMC, 3995 St Marys Rd, Columbus, GA 31907 

What: To the west of St. Mary's Church a thin tree line exists between the church 


and the highway. Test within tree line compared to open. 


Start/end three airbeams in open grass on west side of church. Pick two with closest averages and 
peaks. Leave third at start location. Move two comparable airbeams at similar distances from road 
within trees and the other in open lawn. 

Hypothesis: Thin tree line will not cause difference in particulate levels as compared with open areas. 
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Area P25 wia GA EPD Air Branch: 
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Williams Road 2-28-17 


When: 2/28/2017 15:47-17:00 EDT 

Who: K. Youngquist and Dalton Peters 

Where: Wiliams Road Field across from Shell Gas Station 

What: Cleared field sits beside thick field of trees across the street from gas station and 


off ramp of 1-185. Test within tree line compared to open. 


How: Start/end three airbeams in cleared field near road. Pick two with closest averages 
and peaks. Leave third at start location. Move two comparable atrbeams at similar 
distances from road, within trees and the other in open field. 

Notes: Calm wind might account for highest particulate at control as it was closest to 
street and gas station. Yellow jackets interupted 3rd location test, 3 minutes shorter 
than others. 


Hypothesis: Trees higher particulate matter. trapping gas station and idling car exhaust. 


Particulate Notes 
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Airbeam Location: 
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